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PART  I :  INTRODUCTION 


Overview 


1.  This  report  provides  potential  users  with  a  complete  description  of 
the  method  developed  by  Hanson  and  Kraus  (1985)  for  implementing  the  seawall 
boundary  condition  in  the  shoreline  change  numerical  model.  Example  runs  are 
included  so  that  users  may  test  their  programs.  Computer  programs  written  in 
FORTRAN  77  are  given  and  explained  for  both  explicit  and  implicit  finite- 
difference  numerical  solution  schemes. 

2.  The  governing  principles  for  the  seawall  boundary  condition  are  sum¬ 
marized  in  Part  I.  The  physical  basis  of  the  seawall  boundary  condition  is 
discussed  in  a  general  and  descriptive  way  in  Part  II.  Parts  I  and  II  provide 
background  material  and  can  be  understood  without  knowledge  of  numerical  model¬ 
ing.  Technical  details  of  the  shoreline  numerical  model  and  implementation  of 
the  seawall  boundary  condition  are  given  in  Part  III.  Two  example  calculations 
and  a  discussion  of  numerical  accuracy  and  efficiency  are  given  in  Part  IV. 

The  computer  programs  are  described  in  Part  V  and  listed  in  Appendix  A. 

Purpose  of  Seawalls 

3.  Chronic  erosion  is  found  along  many  portions  of  the  coast  of  the 
United  States  and  other  coasts  of  the  world.  Coastal  erosion  is  caused  by  di¬ 
verse  factors.  These  include  rise  in  mean  sea  level,  increase  in  severity  of 
incident  waves,  change  in  local  magnitude  and  direction  of  incident  waves  (as 
produced,  e.g.,  by  a  newly  installed  coastal  structure),  loss  of  sediment  sup¬ 
ply  from  rivers  and  cliffs,  and  interruption  of  the  local  littoral  drift  by 
structures.  If  the  cause  of  undesirable  erosion  in  an  area  cannot  be  elimi¬ 
nated  or  corrected,  then  buildings,  roads,  and  other  resources  will  eventually 
become  endangered,  and  some  degree  of  shore  protection  must  be  undertaken. 
Chapter  1  of  the  Shore  Protection  Manual  (SPM  1984)  contains  a  detailed  dis¬ 
cussion  of  the  causes  of  coastal  erosion  and  their  remedial  measures. 
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4.  The  shore  can  be  protected  against  erosion  through  the  use  of 
coastal  structures,  nonstructural  procedures,  such  as  beachfill,  or  a  combina¬ 
tion  of  structures  and  nonstructural  methods  (SPM  1984;  US  Army  Corps  of  Engi¬ 
neers  1981).  In  situations  where  extensive  damage  may  occur  because  of  storm 
waves  and  water  intrusion,  or  where  nonstructural  procedures  are  not  feasible, 
then  seawalls,  bulkheads,  and  coastal  dikes  are  commonly  constructed  for  beach 
erosion  control  and  for  preventing  inundation.  If  the  word  "seawall"  is  used 
to  describe  any  man-made  or  natural  object  which  functions  as  a  nonerodible 
barrier  along  the  shoreline,  the  concept  of  "seawall"  encompasses  true  sea¬ 
walls,  coastal  dikes,  storm  surge  barriers,  shore-connected  breakwaters,  bulk¬ 
heads,  revetments,  and  rocky  coastal  cliffs.  A  coast  may  contain  several  such 
seawalls,  and  their  presence  must  be  taken  into  account  when  assessing  the 
long-term  (order  of  years)  evolution  of  the  shoreline. 

5.  It  is  also  necessary  to  estimate  the  impact  of  a  proposed  seawall  in 
the  design  process  for  shore  protection.  Even  a  wide  sandy  beach  cannot  erode 
indefinitely;  at  some  point  in  time  the  beach  material  will  be  exhausted,  and 
permanent  structures  and  resources  will  become  exposed  to  wave  action  and  in¬ 
undation.  In  such  a  situation,  emergency  protective  measures  will  be  taken, 
most  likely  by  the  construction  of  a  revetment,  bulkhead,  or  seawall.  A  nu¬ 
merical  model  of  shoreline  change  must  allow  for  the  real  world  situation  of 
the  ultimate  presence  of  a  seawall. 

Seawalls  and  the  Shoreline  Change  Model 

6.  Numerical  models  provide  a  powerful  means  for  making  quantitative 
estimations  of  shoreline  evolution.  In  particular,  the  so-called  "one-line" 
numerical  model,  originating  from  the  work  of  Pelnard-Considere  (1954),  has 
been  widely  applied  in  recent  years.  Kraus  (in  preparation)  gives  an  anno¬ 
tated  bibliography  of  the  literature  on  one-line  models.  The  term  "one-line" 
typically  refers  to  the  shoreline;  therefore,  this  model  is  often  called  the 
"shoreline"  model.  Despite  the  large  number  of  applications  of  the  shoreline 
model,  representation  of  the  action  of  a  seawall  in  the  model  has  received 
little  attention.  A  seawall  imposes  a  constraint,  or  boundary  condition,  on 
the  solution  (shoreline  position)  obtained  with  the  model. 

7.  The  most  obvious  boundary  condition  imposed  by  a  seawall  is  that 
the  beach  fronting  the  wail  cannot  move  landward  of  it.  Also,  a  seawall 
prevents  the  sediment  contained  behind  it  from  entering  the  littoral  system, 
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thereby  modifying  the  sand  transport  rate  along  the  beach  and  possibly  starv¬ 
ing  the  adjacent  beach  through  the  elimination  of  potential  littoral  material. 
In  an  extreme  case,  if  the  level  of  the  beach  in  front  of  a  seawall  drops, 
waves  will  reflect  from  the  wall  instead  of  dissipating  on  the  beach.  Stand¬ 
ing  waves  can  cause  locax  scour  that  may  temporarily  increase  transport  along¬ 
shore  or  offshore,  until  a  new,  steeper  equilibrium  profile  is  achieved.  The 
integrity  of  the  seawall  may  be  threatened  when  the  beach  elevation  drops. 

8.  In  the  literature,  there  has  been  very  little  discussion  on  repre¬ 
sentation  of  a  seawall  in  the  shoreline  model  or  in  other  models.  Essentially 
all  of  the  work  reported  to  date  has  been  conducted  by  engineers  associated 
with  coastal  engineering  in  Japan.  More  than  25  percent  of  Japan's  34,000-km- 
long  (21,000-mile)  coastline  is  protected  by  seawalls,  coastal  dikes,  armor 
blocks,  and  similar  structures  (Ogawara  1983). 

9.  In  the  early  1970's,  Hashimoto  et  al.  (1971)  discussed  the  behavior 
of  the  longshore  sand  transport  rate  in  front  of  a  seawall  armored  by  blocks. 
They  recommended  the  longshore  transport  rate  be  set  to  zero  if  the  shoreline 
reaches  the  seawall.  Ozasa  and  Brampton  (1980)  treated  the  loss  of  berm  in 
front  of  a  seawall  and  devised  prescriptions  for  introducing  the  action  of  a 
seawall  in  the  shoreline  numerical  model.  In  essence,  their  procedure  also 
consists  of  setting  the  longshore  sand  transport  rate  equal  to  zero  at  calcu¬ 
lation  points  where  the  berm  has  been  removed  and  the  shoreline  has  retreated 
to  the  seawall.  Hanson  and  Kraus  (1980)  gave  a  procedure  in  the  form  of  a 
simple  shoreline  adjustment,  but  this  alone  is  unsatisfactory  because  it  does 
not  conserve  sand  volume.  Tanaka  and  Nadaoka  (1982)  noted  that  the  procedure 
of  setting  the  transport  rate  to  zero  is  not  correct.  They  proposed  two  al¬ 
ternative  methods,  but  unfortunately  their  methods  appear  to  be  arbitrary  and 
incomplete. 

10.  Recently,  Hanson  and  Kraus  (1985)  have  given  an  outline  of  a  well- 
tested  procedure  for  representing  the  action  of  a  seawall  in  balance  with  the 
capability  of  the  shoreline  numerical  model  and  in  accordance  with  three  gen¬ 
eral  principles.  The  present  report  gives  a  complete  description  of  their 
method.  The  physical  reasoning  behind  the  method  is  discussed  in  Part  II. 

The  principles  upon  which  the  method  is  based  are: 

a.  The  shoreline  in  front  of  a  seawall  cannot  recede  landward  of 
the  seawall. 

b.  Sand  volume  must  be  conserved. 
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c.  The  direction  of  sand  transport  alongshore  must  be  preserved  in 
accordance  with  the  natural  direction  of  the  potential  local 
transport. 

11.  Although  the  above-listed  principles  are  easy  to  understand,  their 
implementation  in  a  computer  program  is  considerably  involved,  in  particular, 
for  b  and  c.  The  present  report  describes  well-tested  algorithms  for  imple¬ 
menting  the  seawall  boundary  condition  in  a  general  manner. 

Limitations  of  the  Method 

12.  The  seawall  constraint  should  be  formulated  on  the  same  level  of 
idealization  as  the  shoreline  model.  Thus,  it  is  not  appropriate  in  the  model 
to  consider  wave  reflection  and  sea  bottom  scouring,  and  settling,  flanking, 
and  collapse  of  the  seawall  (for  further  discussion  on  one  or  more  of  these 
topics,  see  Sato,  Tanaka,  and  Irie  1969;  Silvester  1977;  Toyoshima  1979; 

Walton  and  Sensabaugh  1979).  It  should  be  stressed  that  the  procedure  de¬ 
scribed  here  possesses  the  same  limitations  as  well  as  the  same  advantages  as 
the  shoreline  model.  The  seawall  boundary  condition  is  only  valid  to  the  ex¬ 
tent  the  shoreline  model  is  valid. 

13-  One  of  the  most  restrictive  assumptions  made  in  deriving  the  shore¬ 
line  model  is  that  the  beach  profile  remains  unchanged  and  moves  seaward  or 
shoreward  in  parallel  to  itself  (an  assumption  of  equilibrium  of  the  profile). 
In  nature,  however,  if  a  beach  erodes  to  reach  a  vertical  or  nearly  vertical 
seawall,  due  to  wave  reflection  and  scouring,  the  beach  slope  immediately  in 
front  of  a  seawall  is  expected  to  become  steeper  than  the  slope  on  the  adjoin¬ 
ing  beach  without  structures  or  steeper  than  the  original  beach  before  the 
seawall  was  built. 

14.  The  above  discussion  notwithstanding,  examples  can  be  found  in  the 
field  of  the  growth  and  recovery  of  formerly  eroded  beaches  fronting  rough¬ 
faced  sloping  seawalls  (Toyoshima  1979);  nearly  vertical  seawalls  (O'Brien 
1985);  and  even  a  vertical  seawall  (Berrigan  1985a, b).  Because  of  an  apparent 
lack  of  data  in  these  cases,  however,  cause  and  effect  have  not  been  clearly 
distinguished.  That  is,  it  is  not  known  with  certainty  whether  the  seawalls 
promoted  growth  of  the  beaches  in  front  of  them  or  if,  e.g.,  sediment  trans¬ 
port  conditions  changed  to  bring  back  the  beaches  with  no  relation  to  the  sea¬ 
wall,  the  seawall  only  initially  functioning  to  protect  the  land  behind  it.  A 
combination  of  the  two  scenarios  is  also  possible.  Toyoshima  (1979)  states 
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that  a  rough-faced  sloping  and,  ideally,  permeable  seawall  will  promote  recov¬ 
ery  by  dissipating  wave  energy,  similar  to  the  functioning  of  a  natural  beach. 

15.  Based  on  the  results  of  their  laboratory  experiments,  Hattori  and 
Kawamata  (1977)  found  that  a  necessary  condition  for  the  naturally  occurring 
restoration  of  an  eroded  beach  backed  by  a  seawall  is  that  a  surf  zone  exist 
seaward  of  the  wall.  Essentially  the  same  conclusion  had  been  reached  in  an 
earlier  laboratory  study  by  Chestnutt  and  Schiller  (1971).  Clearly,  results 
of  simulations  incorporating  the  seawall  boundary  condition  in  a  shoreline 
model  must  be  interpreted  with  caution. 

16.  In  order  to  account  for  an  alongshore  variation  in  beach  slope,  a 
mechanism  to  allow  for  cross-shore  sand  transport  and  a  more  complicated  nu¬ 
merical  scheme  than  that  used  in  the  shoreline  model  are  required.  Numerical 
models  now  exist  which  account  for  cross-shore  transport  in  a  schematic  way. 
The  "2-line"  model  of  Bakker  (1969)  and  Bakker  et  al.  (1971),  and  the  "N-line" 
model  of  Perlin  and  Dean  (1978,  1983)  are  examples.  Such  models  can,  in  prin¬ 
ciple,  more  realistically  represent  the  beach  slope  in  front  of  a  seawall  than 
can  the  shoreline  model. 

17.  At  present,  however,  these  models,  although  more  sophisticated  than 
the  shoreline  model,  have  limitations  for  engineering  use  stemming  from  lack 
of  knowledge  of  the  physical  mechanism  of  cross-shore  sand  transport.  Numeri¬ 
cal  instability  and  long  computer  run  times  are  the  main  technical  problems 
encountered.  Relatively  short  calculation  time  is  an  appealing  feature  of  the 
shoreline  model.  This  feature,  plus  its  demonstrated  versatility  for  handling 
a  wide  range  of  boundary  conditions,  ensures  the  use  of  the  shoreline  model  as 
an  engineering  tool  in  the  foreseeable  future. 

18.  In  summary  for  this  section,  to  the  extent  that  changes  in  beach 
cross  section  can  be  neglected  in  comparison  to  changes  in  beach  planform,  the 
shoreline  model  is  a  useful  engineering  tool  for  systematically  investigating 
and  estimating  shoreline  evolution  over  time  periods  of  several  months  to  sev¬ 
eral  years.  If  seawalls  are  located  along  the  coast,  because  of  possible  sig¬ 
nificant  changes  in  beach  cross  section,  particular  caution  should  be  exer¬ 
cised  in  interpreting  model  results. 

19.  As  progress  is  made,  it  will  become  desirable  to  incorporate  the 
seawall  boundary  condition  in  models  more  sophisticated  than  the  shoreline 
model.  This  task  may  prove  to  be  difficult.  Experience  and  familiarity  with 
the  implementation  of  the  seawall  boundary  condition  in  the  shoreline  model 
should  provide  useful  guidance. 
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PART  II:  BACKGROUND  FOR  THE  SEAWALL  BOUNDARY  CONDITION 


Action  of  a  Seawall  on  a  Beach 

20.  There  is  remarkably  little  quantitative  information  available  on 
the  behavior  of  real  beaches  backed  by  seawalls.  It  has  been  long  known  that 
under  certain  wave  conditions,  a  vertical  seawall  will  accelerate  erosion  of 
the  beach  in  front  of  it  (see  Russell  and  Inglis  1953;  Sato,  Tanaka,  and  Irie 
1969).  Scour  is  the  primary  cause  of  this  erosion.  Sand  is  scoured  from  the 
sea  bottom  in  front  of  a  vertical  seawall  by  the  standing  wave  system  produced 
by  wave  reflection  at  the  wall.  Any  current,  such  as  the  longshore  current, 
can  then  transport  the  mobilized  sand  out  of  the  area.  If  there  is  a  contin¬ 
ued  net  loss  of  sand  over  a  long  period  of  time,  the  end  result  is  that  the 
beach  in  front  of  the  seawall  can  no  longer  maintain  the  natural  equilibrium 
profile  and  the  beach  slope  will  become  steeper.  Walton  and  Sensabaugh  (1979) 
discuss  this  and  other  processes  believed  to  enhance  erosion  of  beaches  backed 
by  vertical  or  nearly  vertical  seawalls. 

21.  On  the  laboratory  scale,  it  has  been  amply  demonstrated  that  a  sea¬ 
wall  does  not  always  produce  erosion  when  introduced  in  the  active  wave  zone 
of  a  beach  in  equilibrium  with  the  existing  waves.  A  brief  discussion  will 
now  be  given  of  three  experiments  (Dorland  1940,  Chestnutt  and  Schiller  1971, 
and  Hattori  and  Kawamata  1977)  performed  using  sand  beaches  in  two-dimensional 
wave  flumes. 

22.  Dorland  (1940)  used  moderately  steep  waves  in  an  attempt  to  repro¬ 
duce  storm  conditions.  He  placed  a  vertical  seawall  at  the  shoreline  of  a 
beach  which  had  been  allowed  to  attain  equilibrium  under  constant  wave  action, 
scooped  out  part  of  the  bed  in  front  of  the  seawall,  and  then  continued 
applying  the  waves.  In  the  two  such  experiments  performed,  the  outer  bar  moved 
landward  and  the  scooped  out  area  partially  filled  with  sand  from  the  offshore. 
In  a  third  series  of  runs  using  three  sets  of  wave  conditions  varying  cycli¬ 
cally,  Dorland  similarly  found  that  the  scooped  out  beach  partially  recovered. 

23.  Chestnutt  and  Schiller  (1971)  found  that  maximum  erosion  occurred 
if  a  seawall  was  placed  on  an  equilibrium  beach  in  a  "critical"  region  lying 
from  about  0.5  to  0.67  xb  ,  where  xb  is  the  width  of  the  surf  zone,  as 
measured  from  the  shoreline.  When  the  seawall  was  moved  to  a  position  shore¬ 
ward  of  the  critical  region,  the  profile  immediately  seaward  of  the  wall  began 


to  accrete,  i.e.,  the  previously  wave-scoured  region  tended  to  be  filled. 
Chestnutt  and  Schiller  point  out  that  the  surf  zone  width  depends,  in  part, 
on  the  wave  period.  Other  factors  being  the  same,  the  surf  zone  will  be  wider 
for  longer  period  waves.  Therefore,  whether  or  not  a  seawall  will  tend  to 
promote  erosion  or  accretion  depends  on  the  wave  conditions,  which  usually 
have  a  marked  seasonal  variation. 

24.  Hattori  and  Kawamata  (1977)  recorded  beach  profile  changes  on  a 
laboratory  beach  with  and  without  a  vertical  seawall.  For  given  wave  con¬ 
ditions,  the  beach  was  allowed  to  attain  equilibrium  before  introduction  of 
the  seawall.  Incident  wave  steepness  was  varied  for  a  fixed  location  of  the 
seawall  relative  to  the  initial  shoreline.  They  found  the  existence  of  a 
surf  zone  to  be  a  necessary  condition  for  recovery  of  an  eroded  seawall -backed 
beach.  This  result  is  in  agreement  with  the  findings  of  Chestnutt  and  Schil¬ 
ler  (1971).  The  existence  of  a  surf  zone  implies  minimum  wave  reflection  at 
the  seawall.  Hattori  and  Kawamata  also  found  that  the  restoring  wave  con¬ 
ditions  for  a  seawall-backed  beach  are  similar  to  those  for  a  natural  labora¬ 
tory  beach  without  a  seawall. 

25.  Movable  bottom  laboratory  experiments  are  difficult  to  interpret 
because  of  scale  effects,  and  longshore  processes  were  absent  in  the  experi¬ 
ments  under  discussion.  Nevertheless,  a  reasonable  conclusion  to  be  drawn 
from  the  aforementioned  work  is  that  an  eroded  beach  in  front  op  a  seawall 
tends  to  recover  when  the  mean  water  level  is  low,  the  waves  have  mild  steep¬ 
ness,  and  a  sediment  supply  exists  in  the  offshore.  Toyoshima  (1979),  O'Brien 
(1985),  and  Berrigan  ( 1985a, b)  give  examples  of  prototype  beaches  backed  by 
seawalls  which  have  become  stable  or  have  recovered  to  some  degree. 

26.  The  interaction  between  beaches  and  seawalls  is  far  from  under¬ 
stood.  A  focused  and  intensive  field  monitoring  effort  is  definitely  needed 
as  a  first  step  toward  achieving  quantitative  understanding  of  the  influence 
of  a  seawall  on  the  shoreline  and  beach  profile.  Without  data,  quantitative 
understanding  and  numerical  modeling  of  the  processes  involved  will  be  limited 
and  suspect. 


Seawall  at  Oarai  Beach,  Japan 

27.  The  physical  picture  for  the  seawall  boundary  condition  formulated 
by  Hanson  and  Kraus  (1985)  is  based  on  general  observations  of  the  shoreline 
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change  at  the  seawall  located  south  of  Oarai  Harbor,  Ibaraki  Prefecture, 

Japan.  A  location  map  is  given  in  Figure  1.  Shoreline  change  at  this  site 
has  been  extensively  documented  and  numerically  modeled  (Kraus,  Hanson,  and 
Harikai  1985).  There  are  two  seawalls  on  this  north-south  oriented  sandy 
beach  facing  the  Pacific  Ocean.  The  north  seawall  is  a  continuous  massive 
concrete  wall  2  km  (1.24  miles)  long  and  5  m  (16.4  ft)  high  from  base  to 
crown.  Portions  of  the  north  seawall  at  Oarai  are  shown  in  Figures  2  and  3. 
The  face  of  the  north  seawall  is  mildly  curved  outward  and  armor  blocks  have 
sometimes  been  placed  at  the  foot  of  the  wall  when  the  beach  eroded.  The 
south  seawall  is  similarly  constructed  and  800  m  (0.5  miles)  long.  Beach 
change  at  the  north  seawall  has  mainly  been  studied. 

28.  When  the  shoreline  reaches  the  seawall,  the  local  beach  slope 
becomes  slightly  steeper  than  the  typical  nearshore  slope  on  this  coast  (which 
itself  varies  between  approximately  1/50  and  1/70  from  the  beach  face  to  the 
wave  breaker  line).  The  change  in  beach  slope  is  mild  and  appears  to  be  neg¬ 
ligible  for  purposes  of  applying  the  shoreline  model.  No  drastic  alteration 
in  beach  characteristics  occurs  and  the  beach  is  exposed  at  low  tide  (Fig¬ 
ure  3).  At  high  tide,  when  the  shoreline  has  receded  to  the  seawall,  broken 
waves  slap  against  the  face  of  the  wall. 

29.  Although  the  shoreline  may  reach  the  seawall  at  some  location,  it 
has  been  inferred  on  the  basis  of  the  observed  and  modeled  long-term  shoreline 
change  that  sand  moves  alongshore  through  the  area  to  be  deposited  adjacent  to 
a  large  groin  at  Oarai  Harbor  (Kraus,  Harikai,  and  Kubota  1981;  Mizumura  1982; 
Kraus  and  Harikai  1983;  Kraus,  Hanson,  and  Harikai  1985).  Since  alongshore 
variations  in  the  slope  of  the  beach  in  front  of  the  seawall  are  small,  the 
seawall  does  not  appreciably  alter  the  pattern  of  wave  breaking.  A  surf  zone 
usually  exists  in  front  of  the  seawall  and  the  capacity  for  waves  to  move  sand 
alongshore  is  retained.  Sand  is  transported  in  the  direction  of  the  wave- 
induced  longshore  current,  and  the  beach  in  front  of  the  seawall  has  been  ob¬ 
served  to  periodically  erode  and  recover. 

Idealized  Seawall  Boundary  Condition 


30.  From  the  observations  described  above,  Hanson  and  Kraus  (1985)  de¬ 
veloped  the  concept  of  the  idealized  functioning  of  a  seawall  for  use  with  the 
shoreline  model.  They  concluded  that  once  the  shoreline  reaches  a  seawall  at 
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Figure  2.  North  seawall  at  Oarai,  Japan,  May  1980 


Figure  3.  South  end  of  north  seawall  at  Oarai,  Japan,  May  1980 
(Seawall  face  in  lower  right  portion  of  photograph) 


a  particular  location,  sand  cannot  originate  from  that  area.  There  can  be  a 
net  gain,  but  no  net  loss,  for  a  beach  area  in  contact  with  a  seawall  (since 
it  is  assumed  in  the  shoreline  model  that  the  beach  level  does  not  drop  below 
the  water  line  and  that  the  beach  slope  does  not  change).  However,  sand  can 
move  alongshore  through  such  an  area,  passing  into  and  out  of  its  boundaries, 
according  to  the  natural  direction  of  transport.  Sand  can  also  be  deposited 
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in  front  of  a  seawall,  thus  allowing  the  beach  to  recover. 

31.  In  the  shoreline  model,  it  would  be  incorrect  to  set  the  transport 

rate  equal  to  zero  at  a  location  where  the  shoreline  makes  contact  with  a  sea¬ 
wall,  as  done  in  most  previous  treatments.  Rather,  the  transport  rate  should 

be  adjusted  to  allow  calculation  cells  in  contact  with  a  seawall  to  transfer 
sand  in  order  to  conserve  total  sand  volume  and  preserve  the  direction  of  its 
transport. 

32.  On  real  beaches,  sand  is  not  always  transported  in  the  same  direc¬ 
tion  over  the  full  length  of  the  beach.  Changes  in  the  direction  of  transport 

may  be  produced,  for  example,  by  longshore  variations  in  wave  direction  and 
wave  height  as  caused  by  refraction  over  an  irregular  bottom,  or  by  diffrac¬ 
tion  at  structures  and  headlands.  Therefore,  at  one  or  more  areas  along  a 
beach,  it  is  possible  that  a  net  amount  of  sand  is  moving  out  of  the  area. 

The  ways  in  which  this  can  occur,  and  implications  for  shoreline  change  in 
the  presence  of  a  seawall,  are  described  in  the  section  Model  Input  Require¬ 
ments  and  Boundary  Conditions,  in  Part  III. 
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PART  III:  SHORELINE  MODEL  AND  THE  SEAWALL  BOUNDARY  CONDITION 


Shoreline  Model  Review 


33.  The  theory  of  the  shoreline  model  originated  with  Pelnard-Considere 
(1954).  He  assumed  that  the  beach  bottom,  not  necessarily  of  planar  slope, 
always  remains  in  equilibrium  and,  as  a  consequence,  moves  in  parallel  to  it¬ 
self  down  to  a  certain  depth,  herein  called  the  depth  of  closure.  Therefore, 
one  contour,  or  "line,"  is  sufficient  to  describe  changes  in  beach  planform. 
This  line  is  conveniently  taken  as  the  shoreline.  Pelnard-Considere  did  not 
develop  a  numerical  model  but  did  give  closed-form  mathematical  solutions  for 
certain  idealized  cases  and  verified  the  results  through  laboratory  experi¬ 
ments.  Details  of  the  numerical  formulation  of  the  model  may  be  found  in, 
e.g.,  Komar  (1976,  1983),  Le  M£haut£  and  Soldate  (1978)  and  Hanson  and  Kraus 
(1980). 

34.  The  purpose  of  the  shoreline  model  is  to  simulate  long-term  evolu¬ 
tion  of  the  shoreline  or  the  beach  planform.  The  governing  equation  for  the 
shoreline  position  is  obtained  from  the  continuity  equation  for  beach  sediment 
(assumed  to  be  cohesionless  sand).  A  predictive  formula  for  the  sand  trans¬ 
port  rate  is  necessary  to  solve  the  governing  equation.  Sand  transport  and 
the  resultant  shoreline  change  depend  on  the  local  wind,  waves,  and  currents, 
beach  planform,  boundary  conditions,  and  constraints  such  as  the  one  produced 
by  a  seawall.  It  will  be  assumed  here  that  the  longshore  sand  transport  is 
produced  solely  by  obliquely  incident  waves;  other  transport  mechanisms  are 
possible,  such  as  coastal,  tidal,  and  wind-generated  currents. 

35.  In  the  present  work,  it  will  be  sufficient  to  use  the  equation  for 
the  shoreline  position  in  its  most  basic  form: 


ll  liS. 

3t  +  D  3x 


(1) 


where 


t 

D 

Q 


shoreline  position,  m 
time,  s 

depth  of  closure,  m 

volume  rate  of  longshore  sediment  transport,  m^/s 
distance  alongshore,  m 


For  simplicity,  only  longshore  transport  of  sand  is  considered.  It  is 
straightforward  to  generalize  Equation  1  to  formally  include  contributions 
for  cross-shore  transport,  as  well  as  sediment  sources  and  sinks.  An  equation 
given  by  Hallermeier  (1979,  1983)  for  a  limiting  depth  of  sand  motion  in  terms 
of  the  incident  wave  conditions  has  been  recommended  by  Kraus  and  Harikai 
(1983)  for  use  as  the  depth  of  closure  (see  also  Kraus  1984). 

Model  Input  Requirements  and  Boundary  Conditions 

36.  In  order  to  solve  Equation  1,  three  kinds  of  information  are  re¬ 
quired:  (a)  the  initial  location  of  the  shoreline  with  respect  to  some  coor¬ 
dinate  system  (Figure  4)  in  which  the  x-axis  is  oriented  along  the  trend  of 
the  coast  and  the  y-axis  points  offshore,  (b)  an  expression  for  the  longshore 
sand  transport  rate,  Q  ,  and  (c)  boundary  conditions  for  either  y  or  Q  at 
the  two  lateral  ends  of  the  beach.  Of  these,  the  initial  position  of  the 
shoreline  is  readily  obtained  or  assumed. 


37.  The  longshore  transport  rate,  Q  ,  is  usually  calculated  from  the 
"CERC"  formula  (SPM  1984,  Chapter  4): 

o  *  *'  ("*  Cg)b  *•"  2  »bs  (2a) 
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where 

K  =  dimensionless  empirical  coefficient  (of  order  0.4) 

H  =  significant  wave  height,  m 
C  =  wave  group  velocity,  m/s 

O 

9bg  =  angle  of  breaking  waves  to  the  shoreline,  deg 
S  =  ratio  of  sand  density  to  water  density 
a'  =  volume  of  solids/total  volume 

r  =  conversion  factor  from  Root  Mean  Square  (RMS)  to  significant  wave 
height,  if  necessary  (equals  1.416) 

The  subscript  b  indicates  quantities  at  wave  breaking.  The  group  velocity 
at  breaking  is  calculated  from: 


^  ■  & r 

where 

g  =  acceleration  of  gravity,  m/s^ 

t  =  ratio  of  wave  height  to  water  depth  at  breaking,  approximately 
equal  to  0.78 

38.  The  angle  9bg  is  the  angle  of  the  breaking  waves  to  the  shore 
line.  It  is  equal  to  the  difference  between  the  angle  the  breaking  waves 
makes  with  the  x-axis  and  the  angle  the  shoreline  makes  with  the  x-axis: 


where 

9fa  =  angle  of  breaking  waves  to  x-axis,  deg 

39.  Common  lateral  boundary  conditions  are  Q  =  0  at  an  impermeable 
barrier  such  as  a  long  jetty  or  groin,  and  3Q/3x  =0  on  a  beach  that  has  a 
stable  (fixed)  shoreline  position.  The  latter  boundary  condition  on  Q  can 
also  be  expressed  as  3y/3t  =  0  (see  Equation  1). 

40.  In  addition  to  lateral  boundary  conditions,  which  are  necessary  to 
solve  any  problem,  it  is  sometimes  required  to  constrain  the  solution,  i.e., 
restrict  movement  of  the  shoreline  position.  For  example,  the  shoreline  along 
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the  beach  backed  by  a  seawall  cannot  recede  behind  the  wall.  In  this  report, 
the  seawall  constraint  is  referred  to  as  a  boundary  condition  although  it  is 
not  a  boundary  condition  in  a  true  sense. 

41.  Three  terms  will  be  defined  to  distinguish  important  transport  sit¬ 
uations  which  can  occur  at  a  seawall. 

Minus  area  (Figure  5a) 

42.  The  expression  "minus"  area  (minus  calculation  cell  in  the  numeri¬ 
cal  model)  is  applied  if,  at  a  given  time,  sand  is  transported  out  of  both 
sides  of  the  area.  If  a  minus  area  occurs  where  the  shoreline  has  eroded  to  a 
seawall,  then  the  sand  transport  rate  must  be  corrected  in  such  a  manner  as  to 
conserve  sand  volume  and  preserve  direction  of  transport,  in  order  to  pass  in¬ 
formation  about  the  lateral  boundary  conditions.  In  the  method  described  in 
this  report,  transport  rate  corrections  along  the  beach  are  made  in  the  direc¬ 
tion  of  sediment  transport,  i.e.,  in  the  downdrift  direction.  Therefore, 
minus  cells  are  starting  points  for  corrections. 

Plus  area  (Figure  5b) 

43.  If  sand  is  moving  into  an  area  from  both  sides  at  a  given  time, 
this  condition  defines  a  "plus"  area  (plus  calculation  cell  in  the  numerical 
model).  The  terminology  "plus  cell"  describes  the  reverse  situation  of  a 
minus  cell;  consequently,  transport  rate  corrections  end  at  plus  cells  (or  at 
lateral  boundaries). 

Regular  area  (Figure  5c) 

44.  The  most  common  situation  is  for  a  certain  quantity  of  sand  to  en¬ 
ter  one  side  of  an  area  and  for  a  slightly  different  quantity  of  sand  to  leave 
the  area  on  the  opposite  side.  This  is  called  a  "regular"  area  (regular  cell 
in  the  numerical  model).  Sand  volume  and  direction  of  transport  must  be  pre¬ 
served  whether  or  not  there  is  a  local  net  gain  or  net  loss  of  material.  If 
the  shoreline  in  a  regular  area  is  in  contact  with  a  seawall,  no  more  sand  can 
leave  the  cell  than  enters  it.  If  the  converse  occurs,  causing  the  nonphys¬ 
ical  movement  of  the  shoreline  to  a  position  landward  of  the  seawall,  the 
transport  rates  must  be  corrected  in  an  appropriate  manner  to  move  the  shore¬ 
line  position  to  the  seawall. 

45.  If  a  wide  beach  exists  in  front  of  a  seawall,  it  is  not  necessary 
to  distinguish  between  minus  areas,  plus  areas,  and  regular  areas.  These 
three  concepts  become  important  only  when  the  shoreline  makes  contact  with  a 
seawall . 


defined  alternately.  Q-points  define  calculation  cell  faces  and  y-points  lie 
at  the  centers  of  cells.  Subscripts  denote  locations  of  points  along  the 
beach.  Both  Q-grid  points  and  y-grid  points  are  separated  by  a  constant  dis¬ 
tance  Ax  alongshore;  the  distance  between  a  Q-point  and  an  adjacent  y-grid 
point  is  Ax/2  .  Lateral  boundary  conditions  must  be  specified  at  the  ends 
of  the  grid,  e.g.,  at  and  °N+1  •  Alternatively,  it  is  possible  to  spec¬ 

ify  boundary  conditions  at  y^  and  y^  ,  or  impose  a  condition  on  y  at  one 
end  of  the  grid  and  a  condition  on  Q  at  the  other  end. 


47.  For  simplicity,  only  one  seawall  will  be  considered.  Its  beginning 
and  ending  coordinates  on  the  x-axis  are  denoted  by  YSBEG  and  YSEND,  re¬ 
spectively,  as  shown  in  Figure  4.  A  general  y-position  at  the  seawall  is  de¬ 
noted  by  ys^  . 

48.  In  a  standard  explicit  scheme,  Equation  1  is  discretized  as 


y!  =  2B  (Q.  -  Q.  .)  +  y. 

J  l  i  1+1  Ji 


(5) 


where 


B  =  At/(2DAx)  ,  s/m 
At  =  time  step,  s 
Ax  =  space  interval,  m 
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49.  For  notational  convenience,  a  prime  on  a  quantity  will  denote  its 
value  at  the  next  (future)  time  step;  an  unprimed  quantity  is  evaluated  at  the 
present  time  step.  Quantities  at  the  present  time  step  are  known.  In  custom¬ 
ary  notation,  the  next  time  step  is  denoted  by  a  superscript  n  +  1  and  the 
present  time  step  is  denoted  by  a  superscript  n  .  The  customary  notation 
will  be  used  in  certain  applications  to  follow. 

50.  For  the  purpose  of  implementing  a  boundary  condition,  or  con¬ 

straint,  the  explicit  model  is  convenient  since  (a)  only  immediately  neigh¬ 
boring  values  of  and  y^  are  involved,  and  (b)  the  implementation  only 

involves  the  present  shoreline  position  and  present  transport  rates;  no  quan¬ 
tities  at  the  next  time  step  are  used. 

51.  If  the  shoreline  moves  landward  of  the  position  of  the  seawall  at  a 
certain  grid  point,  thus  violating  the  seawall  constraint,  the  longshore  sand 
transport  rate  must  be  corrected  to  conserve  sand  volume.  The  (nonphysical) 
erosion,  or  retreat,  of  the  shoreline  to  a  position  behind  a  seawall,  as  shown 
in  Figure  7,  results  in  a  nonphysical  additional  transport  of  sand  out  of  the 
associated  calculation  cell.  The  transport  rates  at  the  cell  faces  must 
therefore  be  corrected  to  prevent  the  shoreline  from  moving  behind  the  sea¬ 
wall.  The  correction  must  be  made  with  consideration  of  the  direction  of 
transport  at  the  two  faces  of  the  particular  cell  violating  the  seawall  con¬ 
straint.  Only  minus  cells  and  regular  cells  may  require  correction.  The  sea¬ 
wall  constraint  is  never  violated  at  a  plus  cell,  because  the  shoreline  always 
advances  in  a  plus  cell. 

52.  The  calculation  procedure  is  described  in  detail  next.  An  overview 
is  as  follows.  First,  the  transport  rates  along  the  beach  are  calculated  in 
order  to  determine  the  transport  directions  and  to  identify  minus,  plus,  and 
regular  cells.  Then,  as  required,  corrections  start  at  either  a  seawall  bound¬ 
ary  or  the  first  minus  cell  encountered  in  the  search.  After  the  starting 
cell  is  corrected,  corrections  to  regular  cells  are  made  as  necessary  follow¬ 
ing  the  direction(s)  of  the  longshore  transport,  until  either  a  plus  cell  or  a 
lateral  boundary  is  reached.  This  procedure  is  repeated  at  each  time  step. 

Correction  at  a 
minus  cell  (Figure  7a) 

53.  Since  correction  is  necessary,  the  shoreline  position  y!  lies 
behind  the  seawall.  The  general  principle  governing  transport  corrections  is 
that  the  transport  rate  at  a  downdrift  cell  face  should  be  reduced  to  a  value 
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a.  Correction  at  a  minus  cell 


b.  Correction  at  a  regular  cell 

Figure  7.  Conceptual  diagram  showing 
shoreline  and  transport  corrections 
at  minus  and  regular  cells 


that  will  place  the  shoreline  at  the  seawall.  In  a  minus  cell,  the  transport 
rates  at  both  cell  faces  are  directed  outward;  therefore,  both  need  to  be  ad¬ 
justed.  It  does  not  appear  that  the  adjustments  can  be  specified  in  a  unique 
way.  Hanson  and  Kraus  (1985)  calculate  corrected  transport  rates  as  equal 
proportions  of  the  original  rates,  as  follows  (with  corrected  quantities  de¬ 
noted  by  a  superscript  asterisk): 
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(6a) 


*  yi  ‘  ysi 

Q.  =  Q.  — - r1 

i  i  y .  -  y ! 

j  i  j  j_ 


*  y i  -  ySi 

Qi  +  1  =  Qi+1  y.  -  y! 

J  l  j  j_ 


The  logic  behind  Equation  6  is  perhaps  more  clearly  understood  by  rearranging 
terms,  to  give,  for  example, 

Q  Q* 

- ± -  =  - ± -  (7> 

yi  '  yi  yi  "  ySi 

from  which  it  is  seen  that,  whereas  Q.  causes  the  shoreline  to  move  from  y. 

1  * 

to  y|  ,  the  corrected  transport  rate  moves  the  shoreline  from  y.  to 

ysj^  (as  required ) . 

54.  By  substitution  of  Equation  6  into  Equation  5,  with  Qi  and 

*  * 

Qi+^  replaced  by  and  »  respectively,  it  is  verified  that  the 

corrected  shoreline  position  is 


yi  *  ysi 


Since  the  adjustment  was  made  through  use  of  the  continuity  equation,  the  pro¬ 
cedure  conserved  sand  volume. 

Correction  at  a 
regular  cell  (Figure  7b) 

55.  With  corrections  at  the  minus  cell  completed,  adjustments  continue 
for  cells  on  both  sides,  following  the  direction  of  transport.  The  transport 
rate  at  an  updrift  face  will  have  previously  been  corrected  and  should  not  be 
corrected  again.  Assuming  for  the  purpose  of  explanation  that  the  transport 

rate  through  a  particular  cell  is  in  the  positive  x-direction,  the  adjusted 

# 

downdrift  transport  rate  Q.  is  obtained  by  setting  the  new  position  y! 


equal  to  ys^  in  Equation  5,  to  give 


=  y.  +  2B  (q.  -  Q*  .  ) 
i  J  i  \  i  i+I  / 


The  corrected  transport  rate  is  then 


*  ysi  '  yi 

Qi+1  =  Qi  '  ^2B  ~ 
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As  a  simple  check,  insertion  of  Qi+1  from  Equation  10  into  Equation  5  gives 
the  following  desired  result  for  the  corrected  shoreline  position: 

y*  =  ysi  (11) 

Again,  this  is  the  mathematical  statement  of  the  shoreline  constraint. 

*  # 

56.  After  Q  ^  and  are  obtained,  calculation  moves  to  the  next 

grid  point  to  determine  y[+i  in  similar  manner.  Calculation  proceeds  from 
cell  to  cell  in  the  downdrift  direction  along  the  seawall  until  either  a  plus 
cell  or  the  end  of  the  seawall  is  encountered.  If  a  plus  cell  is  encountered 
calculation  of  the  shoreline  position  continues  without  necessity  for  correc¬ 
tion  until  another  minus  point  or  the  end  of  the  seawall  is  encountered. 

57.  On  the  other  side  of  the  original  minus  cell,  where  the  transport 

rate  is  in  the  negative  x-direction,  analogous  corrections  are  made  to  trans¬ 
port  rates  as  described  above,  i.e.,  to  .  This  allows  determination  of 

y|  .j  .  The  calculation  then  proceeds  from  cell  to  cell  in  the  downdrift 
direction. 

58.  The  programs  YSEXP  and  CORRE,  discussed  in  Part  V  and  listed  in  Ap¬ 
pendix  A,  calculate  shoreline  change  with  the  explicit  numerical  scheme  for  a 
beach  backed  by  a  seawall. 


Implicit  Numerical  Model 

59.  Compared  to  the  straightforward  development  for  the  explicit 

scheme,  as  presented  in  the  previous  subsection,  representation  of  the  seawall 
constraint  in  an  implicit  numerical  scheme  is  extraordinarily  complex.  In  an 
implicit  scheme,  values  of  the  new  are  solved  for  simultaneously,  over 

the  whole  grid,  in  terms  of  the  old  and  other  quantities.  Thus,  in 

checking  to  determine  whether  the  seawall  constraint  has  been  violated,  the 
time  level  halfway  between  the  old  and  new  time  levels  is  involved.  In  the 
explicit  method,  transport  rates  of  only  those  cells  in  contact  with  a  seawall 
need  to  be  corrected;  in  the  implicit  scheme,  correction  of  one  cell  will  af¬ 
fect  all  cells  downdrift  (whether  in  front  of  the  seawall  or  not)  and  thus  all 
cells  downdrift  require  correction.  Correction  of  all  downdrift  cells  in¬ 
creases  the  complexity  and  execution  time  of  the  computation. 

60.  As  already  discussed,  the  direction  of  sand  transport  must  be 
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preserved  when  correction  of  the  transport  rate  is  made  to  satisfy  the  seawall 
constraint.  Since,  in  general,  the  transport  direction  can  reverse  along  a 
beach,  in  an  implicit  scheme  the  transport  rate  must  be  solved  for  twice, 
starting  independently  from  each  of  the  two  lateral  boundaries.  This  doubles 
the  number  of  calculations  performed,  even  if  no  corrections  are  required,  and 
greatly  reduces  the  speed  advantage  the  implicit  method  normally  holds  over 
the  explicit  solution  method.  Kraus  and  Harikai  (1983)  discuss  and  compare 
the  relative  efficiencies  of  the  explicit  and  implicit  numerical  schemes  for 
the  shoreline  model  without  inclusion  of  the  seawall  constraint.  A  similar 
comparison  of  relative  efficiency,  including  operation  of  the  seawall  con¬ 
straint,  is  given  in  the  examples  discussed  in  Part  IV. 

61.  The  finite  difference  equations  in  an  implicit  scheme  will  be  de¬ 
rived  for  calculating  shoreline  change  in  the  presence  of  a  seawall.  The  grid 
and  notation  are  the  same  as  those  used  in  the  explicit  scheme,  described  in 
the  previous  subsection.  As  the  starting  point,  Equation  1  is  rewritten  to 
give  equal  weight  to  present  and  future  values: 


_  1  /l  3Q  J_  3Q'  \ 

3t  =  "  2  \D  3x  +  D'  3x  / 


(12) 


In  finite  difference  form,  Equation  12  becomes 


y[  --  B'  (Q’  -  Q!  +  1>  +  yc.  (13) 

where 

ycj  =  yt  +  B  (Qt  -  Qi+1)  (14) 

The  quantity  yc^  can  be  interpreted  as  the  shoreline  position  midway  between 
yi  and  y!^  ;  it  is  known  since  it  only  contains  values  at  the  present  time 
step  and  input  data.  The  quantity  B'  =  At/(2D*  Ax)  differs  from  the  un¬ 
primed  version  in  that  it  contains  the  depth  of  closure  at  the  new  time  step, 
which  can  be  calculated  from  the  new  wave  conditions. 

62.  It  is  possible  to  solve  Equation  13  by  an  iterative  procedure  be¬ 
tween  the  y!^  and  the  Qj  ,  as  done  for  example,  by  Le  Mehaut£  and  Soldate 
(1978).  A  computationally  faster  approach  is  to  express  the  Qj  in  terms 
of  the  yj  through  linearization  of  Equation  2.  Such  a  linearization  is 


expected  to  provide  an  accurate  approximation  under  typical  wave  conditions, 
for  which  the  breaking  wave  angle  is  small  (less  than  30).  The  linearization 
method  was  introduced  by  Perlin  and  Dean  (1978)  for  use  with  the  CERC  formula, 
Equation  2.  The  method  was  extended  by  Kraus  and  Harikai  (1983)  to  account  for 
an  additional  contribution  arising  from  a  systematic  change  in  breaking  wave 
height  alongshore  (Ozasa  and  Brampton  1980),  as  caused,  e.g.,  by  wave  diffrac¬ 
tion.  These  references  should  be  consulted  for  details.  The  final  result  is 
that  the  transport  rate  at  the  new  time  step  can  be  expressed  in  the  form 


Q!  =  E!  (y!  .  ♦  y!)  +  F!  (15) 

l  l  Ji-1  l 

where  Ej  and  Fj  are  functions  of  the  incident  wave  parameters.  Substitu¬ 
tion  of  Equation  13  into  Equation  15  gives  a  tridiagonal  system  of  equations 
for  the  Qj  .  A  tridiagonal  system  can  be  solved  by  an  efficient  standard 
algorithm,  called  the  double-sweep  algorithm.  The  solution  is  based  on  the 
following  recurrence  relation: 


Q!  =  EE!  Q!  ,  +  FF! 

l  l  ^1+1  l 


where 


pr  f  _  _ ± _ 

*1  “  1  +  B!  (2  -  EE! 


i-1 


FF! 

l 


F!  +  E!  (yc.  ,  -  yc.)  +  B!  FF!  . 

l  l  J  1-1  J  l  l  1-1 


♦  B!  (2  -  EE!_ .j ) 


(16) 


(17) 

(18) 


B!  :  B  E! 


(19) 


63.  The  solution  procedure,  prior  to  making  any  corrections  to  account 
for  the  seawall,  is  as  follows: 

a.  Specify  a  boundary  condition  at  i  =  1  in  terms  of  EE!  and 

FF!  .  1 

b.  Solve  Equations  17  and  18  for  i  =  2  to  N,  in  ascending  order. 
This  constitutes  the  first  sweep. 

c.  Specify  a  boundary  condition  for  QjJj+1  . 


d.  Solve  Equation  16  for  i  =  N  to  1  ,  in  descending  order.  This 
step  is  the  second  sweep  through  the  grid. 

e.  Substitute  the  Q|  into  Equation  13  to  obtain  the  new  shore¬ 
line  positions,  y!^  . 

64.  The  shoreline  positions  thus  obtained  at  each  time  step  must  be 
compared  with  the  position  of  the  seawall  to  determine  if  the  seawall  con¬ 
straint  was  violated.  If  so,  then  the  shoreline  position  and  associated 
transport  rates  must  be  corrected.  In  general,  when  making  corrections  to 
satisfy  the  seawall  constraint,  it  is  necessary  to  calculate  the  Q!^  in 
ascending  order,  as  well  as  descending  order,  so  that  transport  corrections 
can  be  made  in  either  direction.  The  above  procedure  must  be  repeated  by 
using  a  recurrence  relation  similar  to  Equation  16,  but  which  allows  calcula¬ 
tion  of  Q|  from  the  boundary  condition  at  i  =  1  .  This  relation  has  the 
form: 

Q!  =  PP!  Q!  ,  -*•  RR!  (20) 

l  l  1-1  l 

The  quantities  PP!  and  RR!  depend  on  PP!  .  and  RR!  .  ,  respectively. 

1  1  1+  I  1+  I 

These  quantities  are  defined  similarly  to  EE|  and  FFj  in  Equations  17  and 
18,  and  will  not  be  written  here.  Expressions  for  these  quantities  and  their 
solution  scheme  can  be  found  in  program  YSIMP,  discussed  in  Part  V  and  listed 
in  Appendix  A. 

65.  The  time  evolution  of  y|  in  the  implicit  scheme  is  shown  pictori- 
ally  in  Figure  8a.  For  comparison,  the  analogous  picture  for  the  explicit 
scheme  is  given  in  Figure  8b.  The  shoreline  positions  y^  are  assumed  to  be 
the  same  in  both  cases.  In  the  implicit  scheme,  it  is  seen  that  both  the  pre¬ 
sent  values  (time  level  n)  and  the  future  values  of  Q  (time  level  n  +  1), 
entering  through  3y/3t  in  Equation  12,  are  used  to  calculate  the  shoreline 
change  from  y^  to  y|  .  Since  the  shoreline  change  rates  3y/3t 

are  constant  during  the  time  increment  At  ,  the  shoreline  change  over  a  time 
step  is  a  straight  line. 

66.  The  shoreline  position  midway  (in  time)  between  y^  and  y!  was 
previously  denoted  as  yc^  .  It  is  seen  that  y-points  lie  on  a  straight  line 
between  two  adjacent  yc-points.  Hence,  yc-points  represent  possible  extremes 
in  shoreline  position.  The  important  implication  of  this  is  that  in  the  im¬ 
plicit  scheme  the  seawall  constraint  must  be  formulated  in  terms  of  the  yc^ 
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since  there  is  a  half  time  step  between  ycj  and  and  a  full  time  step 

between  yc^  and  yc!^  .  In  finite  difference  form,  Equation  21  becomes: 

yc!  =  2B '  (Q!  -  Q'+1)  +  yc.  (22) 

A  major  goal  has  been  achieved  by  arriving  at  Equation  22,  because  the  seawall 
constraint  must  be  formulated  in  terms  of  yc-points.  The  implementation  of 
the  constraint  is  similar  to  that  for  the  explicit  scheme,  and  only  an  outline 
will  be  given. 

Correction  at  a  minus  area 

68.  As  in  the  explicit  scheme,  transport  adjustments  start  at  a  minus 
cell  and  from  there  are  performed  in  the  direction  of  transport.  For  the  mi¬ 
nus  cell  itself,  the  adjustment  resembles  that  expressed  by  Equation  6  and 


reads  as  follows: 


yc,  -  ysi 

Q!  =  2(yc1  -  yl)  °i 


yc,  -  ys{ 

Q*  -  - = - —  n> 

i+I  2(yc.  -  y! )  wi- 


i+1  "  2(yci  -  yj)  i+1 

Substitution  of  these  corrected  values  into  Equation  22,  and  using  Equa¬ 
tion  13,  verifies  that  the  desired  result  has  been  obtained,  i.e., 


(23a) 


(23b) 


yc!  =  yst 


Finally,  the  corresponding  corrected  shoreline  position  is  computed  from  Equa¬ 
tion  13  as 

ys,  +  yc 

r-  \ 


The  corrected  position  is  thus  found  to  lie  halfway  between  the  previous  ex¬ 
tremal  position,  yc^  ,  and  the  seawall. 

Correction  at  a  regular 
cell,  positive  transport 

69.  Corrections  are  made  by  moving  in  the  positive  x-direction.  Since 
the  transport  rate  into  the  cell  has  already  been  adjusted  in  connection  with 
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the  previous  (updrift)  cell,  only  the  transport  rate  out  of  the  cell  must  be 
adjusted  in  order  to  satisfy  Equation  20.  This  equation  contains  information 
about  the  upstream  boundary  condition.  Before  any  adjustments  are  made  at 
cell  face  i  +  1  ,  Equation  20  reads 


Qi+1  =  PPl+1 


Q*  + 
1 


RR!  . 
1+1 


(26) 


where  Q?  is  the  corrected  rate  made  for  the  previous  cell.  This  relation 
holds  unless  the  seawall  constraint  was  violated.  If  so,  then  must  be 

adjusted  by  setting  ycj  equal  to  ys^  in  Equation  22,  thus  giving 

ys.  =  2B '  (Q*  -  Q«+1)  +  ycj  (27) 


This  is  easily  solved  for  the  corrected  transport  rate  for  the  downdrift  cell: 

ys.  -  yc. 

Q?+i  =  -  --V—  <28> 

The  procedure  used  to  arrive  at  Equations  26-28  is  continued  in  the  downdrift 
direction  until  either  a  plus  cell  or  a  boundary  is  encountered. 

Correction  at  a  regular 
cell,  negative  transport 

70.  The  procedure  used  here  for  making  corrections  downdrift,  in  the 
negative-x  direction  (on  the  other  side  of  the  minus  cell),  is  very  similar  to 
the  procedure  described  immediately  above.  The  new  transport  rate  at  cell 
face  i  is  given  by  Equation  16,  i.e., 


Q!  =  EE!  Q*  ,  +  FF!  (29) 

l  l  1+1  l 

Then  the  corrected  transport  is  found  to  be 

ys.  -  yc. 

=  Q!+1  *  ^2B~'  ~  (30) 

This  procedure  is  repeated  downstream  until  a  plus  cell  or  a  boundary  is 
encountered . 
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PART  IV:  EXAMPLE  CALCULATIONS 

General  Comments 

71.  Two  examples  are  presented.  These  hypothetical  situations  demon¬ 
strate  applications  of  the  shoreline  model  with  an  operative  seawall  boundary 
condition  and  allow  checking  of  user  implementations  of  the  programs  given  in 
Appendix  A.  An  attempt  was  made  to  give  semi-plausible  examples  while  also 
preserving  clarity.  This  resulted  in  two  idealized  cases  for  which  most  of 
the  common  structures  and  boundary  conditions  could  be  included.  The  first 
example  is  that  of  an  initially  straight  shoreline  bounded  on  one  side  by  a 
jetty.  The  beach  is  protected  by  the  combination  of  a  detached  breakwater  and 
a  straight  seawall  segment.  The  second  example  is  a  curved  pocket  beach  lying 
between  two  headlands  and  protected  by  a  curved  seawall.  Hanson  and  Kraus 
(1985)  show  results  of  several  other  sample  calculations. 

72.  In  the  examples,  the  wave  field  is  introduced  artificially;  the 
breaking  wave  height  and  breaking  wave  angle  were  fabricated  "by  hand"  to 
achieve  the  desired  trends  in  shoreline  movement  in  order  to  exercise  the  sea¬ 
wall  constraint  algorithms.  The  breaking  wave  data  are  set  in  the  subprogram 
INDATA  which  is  given  in  Appendix  A.  Values  of  the  time  and  space  steps  and 
other  parameters  are  entered  via  FORTRAN  DATA  statements.  The  names  of  param¬ 
eters  and  variables  closely  follow  the  notation  of  the  main  text  of  this  re¬ 
port.  The  important  exceptions  are:  the  angle  "theta,"  denoted  as  "Z,"  and 
the  empirical  coefficient  "K,"  denoted  as  "K1"  in  the  program. 

73.  Both  examples  can  be  run  using  either  the  explicit  or  the  implicit 
numerical  scheme,  programs  YSEXP  and  YSIMP,  respectively,  in  Appendix  A.  In 
the  latter  part  of  these  programs  a  calculation  is  made  to  check  sand  volume 
conservation.  It  can  be  verified  that  volume  is  conserved  to  within  trunca¬ 
tion  error. 

Stability 

74.  Before  proceeding  to  the  examples,  the  stability  properties  of 
the  shoreline  model  are  briefly  reviewed.  It  can  be  shown  (e.g.,  Hanson  and 
Kraus  1 980 ;  Kraus  and  Harikai  1983)  that  for  small  breaking  wave  angles  and 
constant  wave  height,  Equation  1,  together  with  Equation  2,  reduces  to  the 
functional  form  of  the  heat  equation,  the  governing  equation  derived  by 
Pelnard-Considere  (1954).  The  accuracy  and  stability  properties  of  numerical 
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schemes  for  solving  this  equation  are  well  known.  Generally  speaking,  numer¬ 
ical  accuracy  can  be  improved  somewhat  by  taking  a  smaller  time  step  for  a 
given  space  step,  assuming  negligible  numerical  truncation  error.  Increased 
computer  execution  time  is  the  price  paid  for  using  smaller  time  steps. 
Therefore,  one  wants  to  balance  speed  of  the  calculation  with  numerical 
accuracy . 

75.  Numerical  accuracy  should  be  distinguished  from  "physical"  accu¬ 
racy.  Numerical  accuracy  is  a  measure  of  how  well  a  finite  difference  scheme 
reproduces  the  solution  of  a  differential  equation;  physical  accuracy  is  a 
measure  of  how  well  the  differential  equation  (and  the  numerical  solution  if 
one  is  employed)  describes  the  process  of  interest. 

76.  For  an  explicit  scheme,  there  is  a  stringent  limitation  (the 
Courant  condition)  on  the  size  of  the  largest  possible  time  step,  other  vari¬ 
ables  being  held  constant.  For  small  breaking  wave  angles,  in  the  present 
case  this  condition  is 
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The  quantity  Rg  was  called  the  "stability  parameter"  by  Kraus  and  Harikai 
(1983).  Equation  31a  is  an  adequate  indicator  of  stability  in  most  applica¬ 
tions,  since  breaking  wave  angles  are  usually  small.  The  stability  parameter 
gives  an  estimate  of  the  numerical  accuracy  of  the  solution,  with  accuracy 
typically  increasing  for  decreasing  values  of  Rg  . 

Example  1 :  Jetty  and  Detached  Breakwater 


77.  The  initial  condition  is  shown  in  Figure  9a.  The  initially 
straight  2,000-m  stretch  of  beach  is  protected  by  a  shore-parallel,  detached 
breakwater  and  a  seawall  connected  to  a  long  jetty.  The  seawall  is  set  back 
7  m  from  the  initial  shoreline.  The  Jetty  is  assumed  to  be  sufficiently  long 
so  as  to  act  as  a  complete  littoral  barrier.  The  breakwater  is  drawn  in  Fig¬ 
ure  9  to  aid  visual  understanding;  in  actuality,  it  lies  much  farther  off¬ 
shore.  The  seawall  and  breakwater  have  been  constructed  to  prevent  erosion 
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of  the  beach  adjacent  to  the  jetty.  The  beach  on  the  far  left  side  of  the 
figure  is  assumed  to  be  outside  the  area  of  influence  of  the  structures,  and 
therefore  its  position  remains  fixed. 

78.  Waves  arrive  at  the  site  as  shown  in  Figures  9b,  c,  and  d.  In 
these  figures,  the  longshore  distributions  of  the  breaking  wave  height  and 
breaking  wave  angle  are  displayed  in  graphic  form  above  the  related  beach 
planform.  The  local  breaking  wave  height  and  angle  are  mainly  controlled  by 
the  detached  breakwater.  The  shoreline  that  would  result  if  there  were  no 
seawall  is  indicated  by  a  dashed  line. 

79.  Figure  9b  shows  the  result  of  waves  arriving  almost  normal  to  the 
shoreline  for  a  period  of  84  hr.  Convergence  of  waves  behind  the  detached 
breakwater  causes  a  bulge,  or  salient,  to  form.  The  wave  direction  then 
changes,  Figure  9c,  and  waves  arrive  obliquely  from  the  right  for  an  elapsed 
time  of  180  hr.  This  results  in  a  loss  of  sand  on  the  beach  next  to  the 
jetty.  The  seawall  prevents  the  shoreline  from  eroding  farther  landward  imme¬ 
diately  next  to  the  jetty;  the  price  paid  is  that  more  sand  is  removed  from 
along  the  front  of  the  seawall.  Finally,  as  shown  in  Figure  9d,  the  wave  di¬ 
rection  changes  again  and  waves  arrive  obliquely  from  the  left.  The  wave 
shadow  zone  behind  the  detached  breakwater  also  shifts  and  the  potential  re¬ 
gion  for  erosion  moves  to  the  middle  of  the  seawall.  Sand  returns  next  to  the 
jetty,  and  an  eroded  sector  forms  at  the  middle  of  the  seawall. 

80.  Although  differences  in  shoreline  positions  with  and  without  the 
seawall  are  moderate  in  this  example,  by  altering  the  input  wave  conditions 
(e.g.,  by  increasing  the  difference  in  breaking  wave  angle  between  applied 
wave  conditions)  a  much  greater  disparity  in  resultant  shorelines  can  be 
generated. 

Example  2:  Pocket  Beach 

81.  The  initial  shoreline  configuration  for  this  example  is  shown  in 
Figure  10a.  A  curved  pocket  beach  approximately  2  km  long  is  bounded  by  two 
long  headlands  which  contain  the  littoral  transport.  A  curved  seawall  is  lo¬ 
cated  4  m  landward  of  the  initial  shoreline. 

82.  Waves  first  arrive  obliquely  from  the  right  side  of  the  figure  for 
126  hr  to  produce  the  planform  shown  in  Figure  10b.  As  a  result,  beach  mate¬ 
rial  moves  toward  the  left  headland.  The  seawall  has  protected  the  area  on 
the  right  side  of  the  beach,  as  seen  by  the  shoreline  change  that  would  have 
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occurred  without  the  seawall  (dashed  line).  The  incident  waves  then  swing 
in  direction  and  arrive  obliquely  from  the  left  for  138  hr,  as  seen  in  Fig¬ 
ure  10c.  Sand  is  transported  past  the  center  of  the  seawall  to  form  a  wide 
beach  adjacent  to  the  right  headland.  The  beach  planform  in  (c)  is  not  a  mir¬ 
ror  image  of  (b)  because,  although  the  waves  were  mirror  images,  the  initial 
shoreline  conditions  were  different. 

83.  In  Figure  10c,  the  seawall  is  protecting  approximately  half  of  the 
shore,  and  much  of  the  eroded  sector  is  still  located  on  the  right  side.  In¬ 
tuition  might  have  suggested  more  erosion  on  the  leftmost  side  since  the  more 
recent  waves  were  from  the  left.  However,  the  interaction  between  waves  and 
shoreline  is  nonlinear  (Equation  2,  the  sine  dependence),  and  the  calculated 
change  is  different  than  might  be  expected.  Finally,  almost  normally  incident 
waves  arrive  to  the  coast  for  72  hr,  to  give  the  result  shown  in  Figure  lOd. 
The  beach  has  essentially  returned  to  its  initial  planform,  Figure  10a.  A 
beach  again  exists  all  along  the  front  of  the  seawall. 

84.  In  this  example,  the  seawall  protected  the  beach  under  episodes  of 
oblique  wave  incidence,  preventing  excessive  landward  retreat  of  the  shore¬ 
line.  The  seawall  therefore  worked  to  promote  recovery  of  the  beach  (compare 
solid  and  dashed  lines  in  Figure  lOd).  It  should  be  cautioned  that  this  re¬ 
sult  is  partially  an  artifact  of  the  assumption  of  an  equilibrium  (constant) 
profile.  In  nature,  the  beach  profile  in  an  eroded  area  would  probably  become 
steeper  than  the  average  beach  profile;  it  then  might  take  a  longer  duration 
of  the  normally  incident  waves  to  cause  the  beach  to  recover. 

Comparison  of  Accuracy  and  Efficiency  of  the  Explicit 
Scheme  and  the  Implicit  Scheme 

85.  The  configuration  of  Example  2  was  used  to  compare  the  numerical 
accuracy  and  efficiency  of  the  explicit  and  implicit  numerical  solution 
schemes  when  operating  under  the  seawall  constraint.  Although  the  results  are 
necessarily  site-dependent,  experience  has  shown  the  trends  to  be  representa¬ 
tive  and  the  conclusions  qualitatively  correct.  Kraus  and  Harikai  (1983)  gave 
a  similar  comparison  of  explicit  and  implicit  numerical  schemes  for  shoreline 
models  without  the  seawall  constraint. 

86.  The  results  of  the  comparison  are  shown  in  Table  1.  The  wave  input 
used  was  that  in  Figure  10b  and  run  for  120  hr.  The  values  of  key  parameters 
were  the  same  as  in  the  previous  examples:  maximum  wave  height  H  :3m, 


T  =  8  s,  DX  =  50  m,  and  D  =  6  m.  In  the  comparison,  the  time  step,  DT,  was 
varied  and  the  reference  or  standard  case  was  taken  to  be  the  explicit  scheme 
with  DT  =  6  hr .  The  relative  accuracy  with  respect  to  the  reference  result, 

,  where  Ay  is  the  change  in  shoreline  position  between  final  and  initial 
positions,  is  given  at  three  locations  on  the  left  side  of  the  beach. 

Table  1 

Stability  and  Accuracy  of  Explicit  and  Implicit 
Numerical  Schemes  with  an  Operative  Seawall 
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87.  The  results  in  Table  1  are  qualitatively  similar  to  those  given  by 
Kraus  and  Harikai  (1983).  The  explicit  model  is  computationally  faster  than 
the  implicit  model  per  time  step;  however,  larger  time  steps  can  be  taken  with 
the  implicit  model  while  preserving  reasonable  numerical  accuracy,  allowing  a 
potential  overall  speed  advantage.  For  example,  the  implicit  model  with  a 
time  step  of  24  hr  and  stability  parameter  of  2.02  is  about  30  percent  faster 
than  the  reference  explicit  result,  yet  still  has  acceptable  numerical  accu¬ 
racy.  Engineering  judgment  must  be  exercised  on  a  case-by-case  basis  to  de¬ 
cide  if  a  24-hr  time  step  will  give  acceptable  physical  accuracy.  In  a  simi¬ 
lar  comparison  without  a  seawall,  Kraus  and  Harikai  (1983)  found  the  implicit 
model  with  a  6-hr  time  step  to  be  comparable  in  accuracy  and  execution  time  to 
the  reference  explicit  model  with  the  same  time  step.  As  was  discussed  in 
Part  III,  the  implicit  model  suffers  a  loss  in  efficiency  when  the  seawall 
boundary  condition  is  operative. 


PART  V:  EXPLANATION  OF  COMPUTER  PROGRAMS 
General  Comments 

88.  Here,  an  explanation  is  given  of  main  operations  performed  in  four 
of  the  five  FORTRAN  programs  given  in  Appendix  A.  The  programs  are  set  up  to 
compute  the  examples  presented  in  Part  IV.  The  final  shoreline  positions  cal¬ 
culated  in  the  examples  are  given  in  Part  IV  so  that  user  implementations  of 
the  programs  can  be  checked. 

89.  The  programs  constitute  the  foundation  of  a  "1-line  model"  and  cal¬ 
culate  shoreline  change  on  a  beach  backed  by  a  seawall  by  means  of  either  the 
explicit  or  the  implicit  numerical  scheme.  In  order  to  run  the  programs  for  a 
general  case,  wave  information  is  needed  to  calculate  the  longshore  sediment 
transport  along  the  beach  in  question.  Specifically,  the  breaking  wave  height 
and  angle  along  the  beach  are  required.  The  breaking  wave  field  must  be  ob¬ 
tained  from  a  wave  calculation  program  such  as  a  refraction  program  or  from  a 
combined  refraction  and  diffraction  program  if  large  coastal  structures  are 
involved.  It  was  beyond  the  scope  of  this  report  to  include  a  numerical  wave 
model.  The  breaking  wave  field  will  also  be  influenced  by  the  plan  shape  of 
the  beach  (the  so-called  sediment-wave  interaction),  which  changes  with  time. 
Numerical  wave  models  and  their  relation  to  the  shoreline  change  model  are 
discussed  by  Kraus  ( 1 983 ) . 

90.  The  five  subprograms  are  called  by  a  main  program.  Input  wave  data 
for  the  examples  are  fabricated  in  subroutine  INDATA.  The  subroutine  INDATA 
is  elementary  and  will  not  be  discussed.  The  longshore  sand  transport  rate, 
computed  by  means  of  Equation  2,  is  calculated  in  subroutines  YSEXP  (explicit 
solution  scheme)  and  YSIMP  (implicit  solution  scheme).  Shoreline  change  in 
the  presence  of  a  seawall  is  computed  in  subroutines  CORRE  (explicit)  and 
CORRI  (implicit).  These  latter  two  routines  correct  both  the  transport  rate 
and  shoreline  position  as  described  in  Part  III. 

91.  Many  of  the  algorithms  are  repeated  in  the  subroutines.  Comments 
are  given  once  for  each  generic  type  of  algorithm.  For  clarity,  the  programs 
are  arranged  to  calculate  for  only  one  continuous  seawall  of  arbitrary  length 
and  configuration.  They  can  easily  be  generalized  to  handle  any  number  of 
seawalls . 

92.  In  the  explanations,  the  names  of  variables  and  line  numbers  refer 
to  those  in  the  indicated  programs.  Line  numbers  in  parentheses  refer  to  the 
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explicit  program  version.  The  names  of  most  key  variables  in  the  programs  are 
the  same  as  those  used  in  the  main  text  of  this  report.  They  are  again  de¬ 
fined  here  to  make  the  explanation  more  self-contained.  The  programs  them¬ 
selves  contain  a  large  number  of  comment  statements  describing  the  operations 
performed  in  distinct  program  segments. 

Programs  YSEXP  and  YSIMP 

93.  Lines  170-190  (150-170):  These  statements  initialize  basic  pa¬ 
rameters.  YSBEG  and  YSEND  define  the  beginning  and  end  grid  points  of  the 
seawall  (Figure  4),  with  YSBEG  <  YSEND.  The  grid  spacing  is  DX  (in  meters) 
and  the  time  step  is  DT  (in  hours).  NTIMES  specifies  the  number  of  timesteps 
and  IT1  and  IT2  denote  timesteps  when  the  wave  data  are  changed  in  the  exam¬ 
ples.  DENOM  is  the  value  of  physical  quantities  in  the  denominator  of  Equa¬ 
tion  2b,  evaluated  for  quartz  sand.  K1  is  the  empirical  coefficient  (K)  in 
Equation  2b.  The  wave  period  is  denoted  by  T  (seconds). 

94.  Lines  250-310  (240-290):  Specify  initial  shoreline  and  seawall 
positions  for  a  straight  beach  and  seawall. 

95.  Lines  370-450  (350-430):  Specify  initial  shoreline  and  seawall 
positions  for  a  curved  beach  and  seawall. 

96.  Line  570  (530):  Call  in  wave  data  and  renew  as  specified. 

97.  Line  730  (680):  Calculate  closure  depth,  DCL0S,  from  wave 
conditions. 

98.  Lines  790-830  (780-810):  These  lines  specify  boundary  conditions 
for  the  simple  cases  of  a  fixed  beach  position  and  an  impermeable  long  groin 
(jetty,  headland). 

99.  Lines  850-960  (720-750):  Calculation  of  the  longshore  transport 

rate. 

100.  Lines  1080-1150:  In  the  implicit  model,  in  order  to  make  correc¬ 
tions  in  both  directions,  a  reversed  double  sweep  is  necessary.  The  longshore 
transport  rates  in  the  arrays  Q  and  QQ  should  be  equal;  a  checking  proce¬ 
dure  is  provided  to  verify  this. 

101.  Lines  1270  (890):  After  the  shoreline  position  is  calculated, 
each  y^  must  be  checked  to  see  if  it  violates  the  seawall  constraint.  The 
subroutines  C0RRI  and  C0RRE  are  called  to  do  the  check  and  to  correct  the 
shoreline  positions  and  transport  rates  as  necessary. 
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102.  Line  1300  (970):  This  program  segment  is  an  error  checking  calcu- 
lation  to  verify  that  sand  volume  was  conserved.  It  also  accounts  for  sand 
that  may  have  entered  the  system  at  the  boundaries. 


Programs  CORRE  and  CORRI 


103.  Subroutines  CORRE  and  CORRI  are  called  by  YSEXP  and  YSIMP,  respec¬ 
tively.  They  recalculate  the  transport  rate  due  to  the  possible  limited  vol¬ 
ume  of  sand  in  front  of  a  seawall  and  adjust  the  position  of  the  shoreline 
accordingly . 

104.  Line  200  (190):  A  branch  is  made  according  to  whether  the  trans¬ 
port  rate  Q.^  is  less  than,  greater  than,  or  equal  to  zero.  A  branch  is  nec¬ 
essary  because  the  corrections  must  be  performed  in  the  direction  of  sand 
transport. 

105.  Line  (200):  This  and  similar  lines  correspond  to  Equation  5. 

106.  Line  260:  Corresponds  to  Equation  20. 

107.  Lines  270-310  (260-310):  If  the  intermediate  shoreline  position 
YC  (for  the  explicit  scheme,  position  Y)  is  seaward  of  the  seawall,  no  correc¬ 
tion  is  necessary.  If  not,  the  downstream  transport  rate  Q^+i  must  be  cor¬ 
rected  in  order  to  conserve  sand  volume.  The  position  YC  (Y)  is  then  set  to 
the  corresponding  position  of  the  seawall. 

108.  Lines  540-680  (540-680):  Calculate  as  described  above,  but  for 
the  reversed  transport  direction. 

109.  Lines  700-800  (700-800):  Corrections  at  a  minus  point  are  com¬ 
puted.  Sand  cannot  be  generated  in  a  minus  cell  located  at  a  seawall.  There¬ 
fore,  the  transport  rates  at  both  cell  faces  are  corrected  so  that  the  shore¬ 
line  will  not  move  landward  of  the  seawall. 

110.  Lines  820-950  (820-970):  This  program  segment  operates  in  the  same 
manner  as  similar  segments  previously  described,  except  that  here  the  calcula¬ 
tion  is  done  in  order  of  decreasing  index  since  the  transport  is  in  the  nega¬ 
tive  x-direction.  Calculation  starts  at  the  point  to  the  left  (lower  i-values) 
of  the  minus  cell  and  continues  downstream  until  a  plus  cell  is  encountered. 

111.  Lines  970-1100:  After  corrections  are  completed  for  grid  points 
within  the  domain  of  the  seawall,  the  same  procedure  must  be  carried  out  for 
the  unprotected  (unstructured)  parts  of  the  beach,  if  any.  This  step  is  nec¬ 
essary  for  the  implicit  scheme,  since  all  values  of  Q  are  solved  at  once. 

It  is  not  required  in  the  explicit  scheme,  for  which  corrections  are  com¬ 
pletely  determined  point  by  point,  at  the  present  time  step. 
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C#  Program  YSEXP  calculates  shoreline  change  according  to  one  lin* 
C*  theory,  taking  into  account  the  effects  of  a  seawall. 

INTEGER  YSBEG, YSEND 
REAL  K1.KAP1 

DIMENSION  Y(40) , YS(40> , SI 41 > , Z 140) ,H(40) , Y0(40> 

DATA  YSBEG/ 26/ , YSEND/ 40/ , DX/50. /,DT/6. / 

DATA  DENOM/2. 362  / , NTIMES/44/ , N/40/ , IT1 / 1 5/ , I T2/ 31 / 

DATA  Kl/0. 1 2/, T / 8. 0/ ,G/ 9. 806/ .GAMMA/ 0. 78/, RADIUS/ 12000. / 
WRITE!#,*)  ’ #####*#£XPL IC I T  CALCULATION*######’ 

WRITE!#,*)  ’ YSBEG-’ , YSBEG, ’  YSEND- ’, YSEND 

C* 

C# 

C*  Initialize  arrays 
C#  Straight  shoreline 
DO  100  1*1 , N 
9 ! I ) -0. 

Y ( I ) -0 . 

100  CONTINUE 

DO  105  1-1 , N 
YStI)— 7. 

105  CONTINUE 
9  <  N  + 1 ) —0 . 

DCLOS-O. 

GOTO  120 

C#  Curved  shoreline 
DO  110  1-1 , N 

BET -AS IN (FLOAT ( 21  - I ) #DX / RADIUS) 

Y (I) -RADIUS# (1. -COS (BET) ) 

Y0(I)-Y(I) 

110  CONTINUE 

DO  115  I-YSBEG, YSEND 
YS ! I ) — Y ( I ) -4 . 

115  CONTINUE 
120  CONTINUE 


WRITE!#, 10)  ( YS ( I ) , I— 1 , N ) 

KAP1-K1/ ( 16. *DEN0M) 

C*  9L-longshore  transport  rate  over  open  boundary 
C*  9L-0 

DO  200  IT-1, NTIMESM 

IF (IT.E9. 1.0R. IT.E9. IT1.0R. IT.E9. IT2)  IC-1 
C#  Subroutine  INDATA  computes  relevant  input  wave  data 
C*  at  any  desired  time  step. 

IF ( IC . E9. 1 )  CALL  INDATA < IT , IT1 , IT2 , H, Z , N, DT ) 

I F ( I T . E9 . NT IMES+ 1 )  THEN 
IC-2 

I HOURS* (IT-1)#INT(DT) 

WRITE!*, 40) 

WRITE!#,#)  ’FINAL  CONDITIONS  (after  ’ , IHOURS, 
ENDIF 

IF(IC.GE.l)  THEN 
WRITE!*, 40) 

WRITE!*, 30)  (YII),I-1,N) 


hours) ’ 


WRITE!*, 30)  (YII),I-1,N) 

WRITE!#, 40) 

WRITE!*, 20)  !9(I),I»1,N) 

ENDIF 

IF (IC. £8. 2)  GOTO  999 
IC-0 

DCLOS-2. 28*H ( 1 ) -68. 5# (H ( 1 ) /T) *#2/G 
B-DT#3600. / ( 2 . #DCL0S#DX ) 

B2-2. #B 


DO  300  1-2, N 

ZBS-Z ( I ) -ATAN! ( Y ( I ) -Y ( I - 1 ) )  / DX) 

8  !  I ) -H ! I ) *#2*S9RT ( G / GAMMA#H ( I ) ) #KAP I #S IN ( 2#ZBS ) 


730  300  CONTINUE 

760  C#  Boundary  conditions: 

770  C*  Pinned  beach 
780  8(1>=8(2> 

790  C#  Groinis) 

800  C*  8 (1 ) -0. 

810  8<N+1)«0. 

820  C# 

830  IF  (YSBEG.GE.3)  THEN 

840  DO  400  I-l.YSBEG-2 

850  Y(I)-Y(I>-B2*<8<Id) -8(I>  > 

860  400  CONTINUE 

870  END IF 

880  C*  Corraction  of  shoreline  in  -front  of  seawall  if  necessary 
890  CALL  C0RRE ( YSBEG, YSEND, 8, B2, Y, YS ) 

900  IF (YSEND. NE .  N)  THEN 

910  DO  500  I-YSEND+l.N 

920  Y ( I)=Y (I) -B2*<8< Id ) -8 (I) ) 

930  500  CONTINUE 

940  ENDIF 

950  C* 

960  C* 

970  C*  Error  calculation  (DIFF:  closed  boundaries,  AROUT:  open  boundary) 
980  C*  QL-8L+8 ( 1 ) 

990  200  CONTINUE 

1000  999  CONTINUE 

1010  DIFF-0. 

1020  AAREA-0. 

1030  DO  600  1-1, N 

1040  DIFF-DIFF+YOII) -Y(I) 

1050  AAREA-AAREA+ ABS ( YO ( I ) - Y ( I ) ) 

1060  600  CONTINUE 

1070  ERROR-DIFF / AAREA 

1080  C*  AR0UT-8L*0T*3600. /DCLOS-DIFF*DX 

1090  C*  ERROR-AROUT/AAREA 

1100  C * 

1110  C***  Output  *** 

1120  WRITE!*,*) 

1130  WRITE (♦, *)  ’LOST  SAND  VOLUME-’ ,ERR0R*100, ’  X’ 

1140  C*  WRITE!*,*)  ' (8L*DT /D-AREA) / ABSAREA*100- ’ , ERROR* 100, ’ X’ 

1150  10  FORMAT (IX, ’SEAWALL  POSITION’ / (IX, 10F8. 2) ) 

1160  20  FORMATdX,  ’LONGSHORE  TRANSPORT ’/( IX ,  1 0F8,  4)  ) 

1170  30  FORMATdX,  'SHORELINE  POSI TI  ON*  /  ( IX ,  1 0F8. 2)) 

1180  40  FORMAT!//) 

1190  STOP 

1200  END 


SUBROUT I NE  CORRE ( YSBEG , YSEND , 0 , B2 , Y , YS ) 

C*  CORRE  recalculates  transport  rates  (Q>  due  to  limited  sand 
C#  volume  in  Front  oF  a  seawall  and  adjusts  the  shoreline 
C*  position  as  necessary.  Explicit  calculation  scheme. 

C# 

INTEGER  YSBEG, YSEND 
REAL  Q  f  41 ) , Y ( 40 )  ,  YS  ( 40) 


I-YSBEG 

IF ( Q ( I ) .GT.O)  THEN 
Y  ( I  -  1  ) *Y ( I - 1 ) -B2#< 
C# 

C*  Q  positive:  Calc  oF  shor 

C*  as  necessary. 


THEN 

1 )  -B2* (9(1 ) -0(1-1 ) ) 


iith  correction  oF  Q  and 


10  IF(Q(I+1).GE.0>  THEN 

Y  <  I  )*Y  (I )  -B2*(Q(I+1 >  -0(D) 

IF ( Y ( I ) . LT . YS  < I ) )  THEN 
DIFF«YS ( I ) - Y ( I ) 

0(1+1) *0(1*1) -DIFF/B2 
Y(I)-YS(I) 

ENDIF 
1*1  +  1 

IF ( I . EO. YSEND* 1 )  GOTO  100 
GOTO  10 
ENDIF 
K-I 
1*1  +  1 

IF ( I . EQ. YSEND* 1 )  THEN 

Y ( I - 1 >  *Y ( I - 1 ) -B2* (0 ( I ) -0(1-1 ) ) 

GOTO  100 
ENDIF 

IF ( I . EO. YSEND)  THEN 
1*1  +  1 
GOTO  30 
ENDIF 
ELSE 

K* YSBEG- 1 

IF ( YSBEG. EO. 1 )  K*1 
ENDIF 
C# 

C*  Q  negative:  Search  -for  a  minus  point.  IF  absent,  calc  Y 
C#  For  the  right  end  element.  Correct  0  as  necessary. 

C* 

20  IF(G(I+1) .LT.O)  THEN 

I-I  +  l 

IF ( I . EQ. YSEND)  THEN 

IF (0 ( I +1 ) . LE . 0 )  THEN 

Y  ( I )*Y (I ) -B2* (0(1  +  1) -Q ( I )  ) 

IF ( Y ( I ) .LT.YS(I) )  THEN 
DIFF*YS ( I ) -Y ( I ) 

0 ( I >  *0 ( I ) +DIFF/B2 
Y(I)-YS(I) 

ENDIF 
GOTO  30 
ENDIF 
ENDIF 
GOTO  20 
ENDIF 


C*  Minus  point:  Corr  oF  0  out  oF  the  element  iF 
C*  behind  seawall. 

C# 

Yd  )  *Y(I>  -B2# (0(1  +  1 )  -0(D) 

IF ( Y ( I ) .LT.YS(I) )  THEN 


shorel ine  moves 


A4 


750 

760 

770 

780 

790 

800 

810  C* 

820  C*  Ca 
830  C*  po 
840  C* 

850  30 


DIFF«YS ( I ) - Y  < I ) 

QDIFF=8<I+1> -Q(I) 

Q< I >»Q< I) -DIFF/B2*  <  Q ( 1 ) /QDIFF) 
Q<I+1)«Q<I+1) -DIFF/B2* (8(1+1) /QDIFF) 
Y(I)=YS(I> 

ENDIF 


1  c  of  Y  starting  ■from  element  to  the  left  of  minus 
int  or  boundary.  Q  is  negative. 


850  30 

860 
870 
880 
890 
900 
910 

920  40 

930 

940 

950  C* 

960  C*  Ca 
970  C*  po 
980  C* 


DO  40  J-I-l.K, -1 

Y< J)=Y(J) -B2*  (8 ( J ♦ 1 ) -Q ( J ) ) 
IF(Y(J) .LT.YS(J) .AND. J.GE. YSBEG) 
DIFF»YS(J>-Y(J) 

Qi  J) -8 ( J) +DIFF/B2 
Y(J)»YS(J) 

ENDIF 
CONTINUE 
I»I  +  1 

IF ( I . GE. YSEND+ 1 )  GOTO  100 


lc  of  Y  starting  from  element  to  the  right  of  minus 
int  or  boundary.  8  is  positive. 


GOTO  10 
CONTINUE 
RETURN 
END 


/-•.’/  .v.' /.v.y. 


100  C*  Program  YSIMP  is  an  implicit  version  of  program  YSEXP  and 
110  C*  calculates  shoreline  change  according  to  one  line  theory, 

120  C*  taking  into  account  the  effects  o-f  a  seawall. 

130  INTEGER  YSBEG, YSEND 

140  REAL  K1.KAP1 

150  D I HENS I ON  Z < 40) , Y < 40 > , YS < 40 > , Y0 < 40) , Q ( 41 > ,YC0LD(40) ,E(40) ,  F (40) 

160  DIHENSI0N  EP ( 40 ) , FP ( 40 ) , BP ( 40 ) , P ( 41 ) , R (41 ) , QQ ( 41 ) , H ( 40) 

170  DATA  YSBEG/1 /, YSEND/40/ , DX/50. / , DT/6. / 

180  DATA  DENOM/2. 362  / , NTIHES/56/ , N/40/ , IT1 /22/ , IT2/45/ 

190  DATA  Kl/0. 12/, T/8.0/,6/9. 806/, GAMMA/0. 78/, RADIUS/ 12000. / 

200  WRITE!*,*)  ’ ******* I HPL I C I T  CALCULATION*******’ 

210  WRITE!*,*)  ’YSBEG*’ , YSBEG, ’  YSEND* ’, YSEND 

220  C* 

230  C*  Initialize  arrays 
240  C*  Straight  shoreline 
250  DO  100  I-l.N 

260  Q(I)»0. 

270  Y  1 1 ) =0 , 

280  100  CONTINUE 

290  DO  105  I-YSBEG, YSEND 

300  YS(I>*-7. 

310  105  CONTINUE 

320  Q(N+l)-0. 

330  D0LD=0. 

340  DCLOS-O. 

350  C*  GOTO  120 

360  C*  Curved  shoreline 
370  DO  110  1*1,  N 

380  BET=ASIN (FLOAT ( 21 - I ) *DX/RADIUS) 

390  Y ( I ) =RADIUS* (l.-COS(BET)) 

400  YO  < I ) — Y ( I ) 

410  110  CONTINUE 

420  DO  115  I-YSBEG, YSEND 

430  YS ( I ) *Y ( I ) -4 . 

440  115  CONTINUE 

450  120  CONTINUE 

460  C* 

470  WR I TE ( *, 10)  (YS ( I ) , 1*1 , N) 

480  KAP 1 *K 1/(16. *DEN0M ) 

490  C*  QL*longshore  transport  rate  over  open  boundary 
500  C*  QL*0. 

510  C*  C*correction  term  in  continuity  calculation 

520  C*  C-1.0 

530  DO  200  IT-1, NTIMES+1 

540  IFIIT.EQ. 1.0R. IT.EQ. 1T1.0R. IT.EQ. IT2)  IC=1 

550  C*  Subroutine  INDATA  computes  relevant  input  wave  data 
560  C*  at  any  desired  time  step. 

570  IF(IC.EQ.l)  CALL  INDATA ( IT , IT1 , IT2 , H, Z , N, DT ) 

580  IFIIT.EQ. NTIMES+l)  THEN 

590  IC-2 

600  I HOURS* ( IT- 1 ) #INT ( DT) 

610  WRITE!*, 40) 

620  WRITE!*,*)  ’FINAL  CONDITIONS  (after  ’ , IHOURS, ’  hours)’ 

630  ENDIF 

640  IF(IC.GE.l)  THEN 

650  WRITE!*, 40) 

660  WRITE!*, 30)  (Y(I),I-1,N) 

670  WRITE!*, 40) 

680  WRITE (*,20)  (Q(I),I-1,N> 

690  ENDIF 

700  IF ( IC. EQ. 2)  GOTO  999 

710  IC-0 

720  DOLD-DCLOS 

730  DC LOS* 2 . 28*H ( 1 ) -68.5*<H( 1 ) /T)**2/G 

740  B»DT*3600. / ( 2. *DCL0S*DX ) 
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rw  irw 7*  v™  "v -  ■%. 


'-  ’1 


730 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 

1130 

1140 

1150 

1160 

1170 

1180 

1190 

1200 

1210 

1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 


B0LD-DT*3600. /  <2.*D0LD*DX> 

YCOLD ( 1 ) =Y ( 1 ) +B0LD*  < Q ( 1 ) -8(2) ) 

C*  Boundary  conditions: 1 
C*  Groin  causing  8<1)»0. 

E(l>-0. 

F ( 1 >  =0. 

C*  Finned  beach  as  8(11=8(2) 

C*  E ( 1 ) =1  . 

C*  F ( 1 ) =0 . 

C* 

DO  300  1-2, N 

YCOLD  d  > -Y  < I > +B0LD# (8  (I > -8 ( I  +  l >  ) 

ZS=ATAN < (Y ( I) -Y (1-1 ) ) /DX> 

Z2«2.*Z(I> 

PWR-H<I>**2*S8RT<G/GAMMA*H(I) ) 

EP  < I ) »PWR*KAP1*2*C0S ( Z2) *<C0S ( ZS) > **2/DX 
FP(I)=PWR*KAP1*SIN(22!*<2*(C0S<ZS>  >**2-1. ) 

BP ( 1 > -B*EP ( I ) 

DEN=1 . +  BP ( I ) * < 2. -E (I -1 )  ) 

E(I) -BP ( I ) /DEN 

F (I)-(FP(I) +EP ( I ) * (YCOLD ( I -1 ) -YCOLD ( I ) ) +  BP (I ) *F ( I - 1 ) ) /DEN 
300  CONTINUE 

C#  Boundary  condition  2:  groin 
8 ( N+ 1 ) —0 . 

C* 

DO  400  I -N ,1,-1 

Q(I)»E(I)*8(I+1)+F(I) 

IF(IT.EQ.l)  THEN 

YCOLD ( I > -Y ( I ) +B* (8 ( I ) -8 ( I ♦ 1 )  ) 

ENDIF 

400  CONTINUE 

C***  Reversed  double  sweep  *** 

C*  Boundary  conditions  3:  groin 
P<N+1)»0. 

R (N  +  l ) =0. 

C* 

DO  500  I-N,2, -1 

P  ( I )  -BP  ( I )  /  (l.iBP(II*(2.'P(Mm 

R< I)-(FP(I) +  EP(I )*( YCOLD <1-1 > -YCOLD ( I ) ) +BP ( I >  *R ( I ♦ 1 ) ) / 

L  <1.*BP(I)*(2.-Pd  +  1>>> 

500  CONTINUE 

C*  Boundary  condition  4  (alt  1:  closed  boundary,  alt  2:  open) 

88(1 )-0. 

C#  8Q(1 )-R(2) / (1 . -P<2) ) 

C* 

DO  530  1-2, N»1 

88 ( I ) -P ( I >  *88 ( I - 1 ) +R ( I ) 

CHECK-ABS (88 ( I ) -8 ( I ) ) 

C*  IFfCHECK.GT. 0.0003)  WRITE!*,*)  ’TRANSPORT  CALC.  DIFFER’ 

550  CONTINUE 

C* 

C*Correction  of  shoreline  in  -front  o-f  seawall  if  necessary 
CALL  CORRI < YSBEG , YSEND ,8,B,YC0LD,E,F,P,R,Y,YS,N) 

C* 

C* 

C*  Error  calculation  (DIFF:  closed  boundaries,  AROUT :  open  boundary) 
C*  IF(IT.EQ.NTIMES)  C-0.5 

C*  8L»8L«-C*8  <  1 ) 

200  CONTINUE 
999  CONTINUE 
DIFF-0. 

A AREA-0. 

DO  600  I-  1 , N 

DIFF»DIFF+Y0(I) - Y ( I ) 

AAREA-AAREA+ABS ( YO ( I ) -Y  (  I ) ) 

600  CONTINUE 
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1410  ERR0R=DIFF /A ARE A 

1420  C*  AR0UT=9L*DT#3600. /DCL0S-DIFF*DX 

1430  C*  ERR0R»AR0UT /AAREA 

1440  C* 

1450  C***  Output  *** 

1460  WRITE  <*, *) 

1470  WRITES*,*)  ’LOST  SAND  VOLUME5*’ ,ERR0R*100, ’  %’ 

1480  C*  WRITE!*,*)  ’ (QL*DT/DCLOS) /ABS (AREA)  *100= ’ , ERR0R*1 00 
1490  10  FORMAT ( IX,  ’SEAWALL  POSI T I ON ’ / ( 1 X , 1 0F8 . 2 )  ) 

1500  20  FORMAT  I IX,  ’LONGSHORE  TRANSPORT ’/ 1 1 X, 10F8 . 4 ) ) 

1510  30  FORMAT ( IX, ’SHORELINE  POSITION’ /( IX,  10F8. 2) ) 

1520  40  FORMAT!//) 

1530  STOP 


"~r-~3rT  'v.  »  - 


SUBROUTINE  CORR I ( YSBEG , YSEND ,9,B,YC0LD,E,F,P,R,Y,YS,N) 
C*  CORRI  recalculates  transport  rates  <9)  due  to  limited  sand 
C*  volume  in  front  of  a  seawall  and  adjusts  the  shoreline 
C*  position  as  necessary.  Implicit  calculation  scheme. 

C* 

INTEGER  YSBEG, YSEND 

REAL  9(41) , Y ( 40 ) ,YS(40) ,YC0LD(40) 

REAL  E<40) ,F<40> ,P(41) ,R(41 ) 

C* 

I-YSBEG 

IF  I Q ( I ) . GT . 0)  THEN 


C*  9  positive*.  Calc  of 
C#  as  necessary. 


shoreline  with  correction  of  9  and 


10  IF (9(1*1 ) . GE. 0)  THEN 

9(1*1) »P (1*1) *9 (I) *R (1*1) 

YC=2*B*(9(I) -9(1*1)  )*YC0LDd) 

IFIYC.LT. YS(I) )  THEN 

D9= ( YS ( I ) -YCOLD ( I ) )/ (2*B> 

9  ( I *1 ) =9 ( I ) -D9 
ENDIF 

Y (I ) =B#  (9 < I ) -9(1*1)) +YC0LD  < I ) 

1  =  1*1 

IF ( I . E9 . YSEND* 1 )  GOTO  100 
GOTO  10 
ENDIF 
K=  I 
1  =  1*1 

IFII.E9. YSEND+1)  THEN 

Y (1-1 )»B* (9(1-1 ) -9 ( I ) ) *YC0LD(I-1 ) 

GOTO  100 
ENDIF 

IF (I. EG. YSEND)  THEN 
1  =  1  +  1 
GOTO  30 
ENDIF 
ELSE 

K-YSBEG 

C* 

C*  9  negative:  Search  for  minus  piont.  If  absent,  calc  Y 
C#  for  right  end  element.  Correct  9  as  necessary. 

C* 

ENDIF 

20  IF (9 ( I *1 ) . LT. 0)  THEN 
1  =  1*1 

IF (I. EG. YSEND)  THEN 

IF ( 9 ( I +1 ) . LE. 0)  THEN 

YC=2*B*<Q<I > -9(1*1)) +YC0LD ( I ) 

IFIYC.LT. YS(I) )  THEN 

DG= ( YS ( I ) -YCOLD ( I ) ) / (2*B> 

Q ( I ) =9 ( I ♦ 1 ) +D9 
ENDIF 

Y  < I ) »B*(9(I) -9(1*1)) +YC0LD  1 1 ) 

GOTO  30 
ENDIF 
ENDIF 
GOTO  20 
ENDIF 
C# 

C*  Minus  point*.  Corr  of  9  out  of  the  element  if  shoreline  moves 
C*  behind  seawall. 

C* 

YC»2*B#  (9(1)  -9(I*D)  +  YCOLD  ( I ) 

IFIYC.LT. YS(I) )  THEN 


750  D8«(YS(  I) -YCOLDd  )  )  /  <2*B> 

760  8DIFF*8(I)  -JIM) 

770  8<I>-a<I)*D8/8DIFF 

780  8(1  +  1)  ■Qd  +  l) *D0/8DIFF 

790  ENDIF 

800  Y ( I ) -B* ( 8 (I > -8 (I ♦ 1 ) > ♦ YCOLD ( I ) 

810  C* 

820  C*  Calc  o-f  Y  starting  -from  element  to  the  left  o+  minus 
830  C#  point  or  boundary.  8  is  negative. 

840  C* 

850  30  DO  40  J»I-1,K,-1 

860  8(J>«E(J>*8(J+1> +F(3) 

870  YC-2*B*(8<J)-8<J+1) >+YC0LD(J) 

880  1F(YC.LT. YS(J) .AND. J.GE. YSBEG)  THEN 

890  DQ« (YS< J) -YC0LD< 3) > / <2*B) 

900  8(J)«8(J+1)+D9 

910  ENDIF 

920  Y< J)«B#<8< J) -Q( J  +  l ) ) +YC0LD  <  J  > 

930  40  CONTINUE 

940  1*1+1 

950  IFd.GE. YSEND+1)  GOTO  100 

960  C* 

970  C#  Calc  o-f  Y  starting  ■from  element  to  the  right  o i  minus 
980  Cft  point  or  boundary.  9  is  positive. 


990  C* 
1000 

GOTO  10 

1010 

100 

CONTINUE 

1020 

DO  110  I-YSBEG-1 ,1,-1 

1030 

8  (I >  =E (I )  #9  <1 ♦ 1 ) +F  <1 > 

1040 

Y  < I ) *B*  <Q<I)-Q<I  +  1> ) +YC0LD  < I ) 

1050 

110 

CONTINUE 

1060 

IF(YSEND.NE.N)  THEN 

1070 

DO  120  I«YSEND+1,N 

1080 

0<I  +  l)*Pd  +  l> *8 (I ) +R (I ♦ 1 > 

1090 

Y ( I ) -B* (8(1) -8(1+1) ) ♦ YCOLD ( I ) 

1100 

120 

CONTINUE 

1110 

ENDIF 

1120 

RETURN 

1130 

END 

100  SUBROUTINE  INDATA ( IT , IT1 , IT2, H, Z, N, DT ) 

110  C*  SPECIFIES  WAVE  HEIGHTS  AND  ANGLES  AT  SPECIFIED  TIME  STEPS 
120  DIMENSION  H(40) ,2(40) ,A(40) 

130  C* 

140  C*IDUM-0: 

ISO  C*Detached  breakwater  version.  The  program  gives  representat i ve 
160  C-Hwave  data  <H,Z)  simulating  effect  of  shore  parallel  1  detached 
170  C*br eakwater  16»DX  offshore  and  running  from  1*12  to  1  =  28. 

180  Cftlnital  beach  is  straight  line. 

190  C* 

200  C#IDUM*1: 

210  C*Represents  an  initially  circular  beach  with  no  offshore  structures. 
220  C* 

230  IDUM-1 

240  DTR-3. 141393/ 180. 

250  IF(IT.EQ.ITl)  GOTO  20 

260  IF ( IT . EQ, IT2)  GOTO  30 

270  C* 

280  C*UAVE  ANGLES 
290  C **#**#•*■*#•**# 

300  C-HCase  1:  Unaffected  breaking  angle  =  0  deg. 

310  IFIIDUM.E8.1)  GOTO  14 

320  DO  10  1=1,12 

330  Z(I)*0. 

340  10  CONTINUE 

350  DO  11  1=13,19 

360  Z(I)=Z(I-l)+10./7. 

370  11  CONTINUE 

380  Z(20)=-10. 

390  DO  12  1=21,28 

400  Z(I)=Z(I-l)+10./8. 

410  12  CONTINUE 

420  DO  13  1=29, N 

430  Z(I)»0. 

440  13  CONTINUE 

450  GOTO  50 

460  C*A1  ter  nat  i  ve  case  1:  Una-f-fected  angle  =  -20  deg. 

470  14  CONTINUE 

480  DO  15  1=1, N 

490  Z  <  I )  *-20.  IfFLOAT  (40  - 1 )  /40. 

500  15  CONTINUE 

510  GOTO  50 

520  C* 

530  C#Case  2:  Unaffected  breaking  angle  =  -10  deg. 

540  20  CONTINUE 

530  IFIIDUM.E8. 1)  GOTO  27 

560  DO  21  1=1,8 

570  Z(I)*-10. 

580  21  CONTINUE 

590  DO  22  1=9,14 

600  Z(I)*Z(I-l)+10./6. 

610  22  CONTINUE 

620  DO  23  1=15, 18 

630  Z ( I ) *Z ( I - 1 ) +2. 5 

640  23  CONTINUE 

650  DO  24  1=19,20 

660  Z ( I ) =Z ( I -1 > -3. 

670  24  CONTINUE 

680  DO  25  1*21,23 

690  Z  < I >  *Z ( I -1 ) -2. 5 

700  25  CONTINUE 

710  DO  26  1*24,40 

720  Z  < I » •- 10. 

730  26  CONTINUE 


All 


740  GOTO  60 

750  C4A1 ternat i ve  cast  2:  Unaffected  angle  *  20  deg 

760  27  CONTINUE 

770  DO  28  1=1,  N 

780  A(I)=-2<Nd-I)/DTR 

790  28  CONTINUE 

800  DO  29  1-1,  N 

810  2 ( I ) =A  < I ) 

820  29  CONTINUE 

830  GOTO  60 

840  C4 

350  C4Ca se  3:  Unaffected  breaking  angle  =  15  deg. 

860  30  CONTINUE 

870  IF(IDUM.ES.l)  GOTO  35 

880  DO  31  1*1, 18 

890  2 ( I ) =15. 

900  31  CONTINUE 

910  DO  32  1=19,28 

920  2(1) “2(1-11-1. 5 

930  32  CONTINUE 

940  DO  33  1=29,33 

950  2 ( I ) =2 (I -1 ) *2. 5 

960  33  CONTINUE 

970  DO  34  1=34,40 

980  2(11=15. 

990  34  CONTINUE 

1000  GOTO  70 

1010  C4A1 ternat i ve  case  3:  Unaffected  angle  =  0  deg 
1020  35  CONTINUE 

1030  DO  36  1*1, N 

1040  2(11=0. 

1050  36  CONTINUE 

1060  GOTO  70 

1070  C4 

1080  C4WAVE  HEIGHTS 

1090  C  4444444444444 

1100  C4Case  1 

1110  50  CONTINUE 

1120  IFdDUM.EQ.il  GOTO  55 

1130  DO  51  1=1,10 

1140  H ( I ) =1 . 50 

1150  51  CONTINUE 

1160  DO  52  1=11,20 

1170  H ( 1 1 =H ( I - 1 1 -0. 1 

1180  52  CONTINUE 

1190  DO  53  1=21,33 

1200  H ( 1 1 «H ( I- 1 1 ♦ 1 . 35/14. 

1210  53  CONTINUE 

1220  DO  54  1=34,40 

1230  H ( 1 1 =1 . 85 

1240  54  CONTINUE 

1250  GOTO  100 

1260  C4A 1  ter nat i ve  case  1 

1270  55  CONTINUE 

1280  DO  56  1*1, N 

1290  H ( I ) =3. 0 -FLOAT (1-1)4.5/40. 

1300  C4  H ( I ) =2. 5 

1310  56  CONTINUE 

1320  GOTO  100 

1330  C4 

1340  C 4-Case  2 

1350  60  CONTINUE 

1360  IFdDUM.EQ.  1)  GOTO  65 

1370  DO  61  1=1,2 

1380  H ( I ) =1 . 85 

1390  61  CONTINUE 


A  1  2 


1400  DO  62  1-3, 16 

1410  Hd)-Hd-1>-1. 33/14. 

1420  62  CONTINUE 

1430  DO  63  1-17,29 

1440  H(I)-Hd-l) +1.33/14. 

1430  63  CONTINUE 

1460  DO  64  1-30,40 

1470  H ( I ) -1 . 85 

1480  64  CONTINUE 

1490  GOTO  100 

1500  C#A1 ter nat i ve  case  2 
1510  63  CONTINUE 

1520  DO  66  1-1 , N 

1330  A  < I ) »H  <N+ 1 - 1 ) 

1540  66  CONTINUE 

1350  DO  67  I— 1 , N 

1560  H ( I ) —A ( I ) 

1570  67  CONTINUE 

1580  GOTO  100 

1390  C* 

1600  C#Case  3 

1610  70  CONTINUE 

1620  IFdDUM.EQ.  1)  GOTO  73 

1630  DO  71  1=1, 12 

1640  HID-1.0 

1650  71  CONTINUE 

1660  DO  72  1-13,24 

1670  Hdl-H(I-l)-0.3/12. 

1680  72  CONTINUE 

1690  DO  73  1-25,37 

1700  Hd)-H(I-l)+l.  33/14. 

1710  73  CONTINUE 

1720  DO  74  1-38, N 

1730  H ( I ) — 1 . 85 

1740  74  CONTINUE 

1730  GOTO  100 

1760  C*A 1  ter nat i ve  case  3 
1770  75  CONTINUE 

1780  DO  76  1-1,20 

1790  A ( I ) -H ( 2*1 ) 

1800  76  CONTINUE 

1810  DO  77  1-21, N 

1820  A ( I ) —A ( N- I ♦ 1 ) 

1830  77  CONTINUE 

1840  DO  78  I — 1 , N 

1850  H ( I >  —A  < I ) 

1860  78  CONTINUE 

1870  C* 

1880  C#OUTPUT 
1890  C* 

1900  100  CONTINUE 

1910  I HOURS- ( IT- 1  HUNT (DT) 

1920  WRITE (4,600) 

1930  IF ( IT. E9. 1 )  WRITE)#,#)  ’ INITIAL  CONDITIONS’ 

1940  IFdT.EQ.  IT1.0R.  IT.EQ.  IT2)  THEN 

1930  WRITE ( # , # )  ’CONDITIONS  AFTER’ , IHOURS, ’HOURS’ 

1960  ENDIF 

1970  WRITE)#,#)  ’WAVE  HEIGHTS’ 

1980  WRITE)#, 602)  (H(I),I-1,N) 

1990  WRITE)#, 600) 

2000  WRITE)#,#)  ’WAVE  ANGLES’ 

2010  WRITE)*, 602)  (2d),I-i,N) 

2020  WRITE)#, 600) 

2030  DO  200  1-1 ,  N 

2040  2 (I ) -Z (I ) *DTR 

2050  200  CONTINUE 

2060  600  FORMAT)//) 

2070  602  FORMAT  (dX,10F8.2)) 

2080  RETURN 

2090  END 

END  OF  FILE 
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Volume  of  solids/total  volume 
Subscript  denoting  breaking  condition 
Wave  group  velocity,  m/s 
Depth  of  profile  closure,  m 

p 

Acceleration  resulting  from  gravity,  m/s 
Wave  height,  m 

Subscript  denoting  position  alongshore 

Dimensionless  empirical  coefficient  in  the  longshore  sediment 
transport  rate  formula 

Superscript  denoting  time  level 

Total  number  of  calculation  cells  in  the  model 

Total  volumetric  longshore  sediment  transport  rate,  m-Vs 

Conversion  factor  from  RMS  to  significant  wave  height 

p 

Stability  parameter,  nr/s 

Ratio  of  density  of  solids  to  density  of  water 
Time,  s 

Position  alongshore,  m 

Position  on-offshore;  shoreline  position,  m 

Extremal,  internally  calculated  shoreline  position,  used  in  the 
implicit  numerical  solution  scheme,  m 

Position  of  seawall  on-offshore,  m 

Ratio  of  wave  height  to  water  depth  at  breaking 

Time  step,  s 

Space  step  alongshore,  m 

Space  increment  on-offshore,  m 

Angle  of  breaking  waves  to  the  shoreline,  deg 

Angle  of  breaking  waves  to  the  x-axis,  deg 

Superscript  denoting  a  corrected  value 

Superscript  denoting  a  quantity  at  next  time  step 


